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In the application of the thermo-electric method of measure- 
ment to the study of the light-curves of variable stars we are con- 
cerned chiefly with relative measurements, and difficulties incident 
to the determination of the ‘‘absolute scale” are of less importance 
in this field. Since in a given series of measurements the compari- 
son stars and variable will be compared image for image, the uncer- 
tain correction to plate center will be largely eliminated. Again, 
we should expect the larger errors due to progressive change in 
atmospheric transparency to be minimized, since presumably 
variable and comparison stars would be similarly affected image 
for image. 

The eclipsing variable U Cephei was selected for measurement 
by the new method for several reasons: the general character of 
the light-curve is well known, elements have been determined, and 
the period is short enough to enable a large portion of the curve to 
be obtained from a single plate. 

Variable-star observers will recall that the light-curve of 
U Cephei has been determined by numerous investigators with 

325 











326 HARLAN TRUE STETSON 


somewhat conflicting results, indicating the possible existence of 
peculiarities yet to be determined. As early as 1884 observations 
by Wilsing' and also by Knott? indicated that the minima observed 
in the spring of the year appeared lower than those deduced from 
autumnal observations by two or three tenths of a magnitude. 
The period of the star being nearly 2} days (more exactly 24493), 
it is obvious that in lower latitudes two successive minima may 
never be seen from the same station. A given observer must 
content himself with observing aliernate minima. Should we regard 
the period as double the foregoing value, the autumnal minima 
would then occupy parts of the light-curve opposite in phase to 
those corresponding to the minima observable in the spring of 
the year. In 1889 Chandler’ published a light-curve showing the 
same peculiarity, indicating also a more rapid rate of increase than 
of decrease, corroborating the results of earlier observations by 
Schmidt. In Chandler’s curve the lowest point is reached first, and 
a short inflection follows near the median of the curve. The reality 
of the seasonal variation, as well as that of the minor fluctuations, 
is regarded by Chandler as open to question. Bohlin,‘ working 
at Upsala from observations made in 1896, secured no foundation 
for any seasonal discrepancies, but obtained a light-curve with a 
slight elevation of about one-tenth of a magnitude at mid-minimum. 

From an extended study of a large number of observations, 
Yendell’ in 1903 came to the conclusion that all such irregularities 
were illusory and attributable to physiological causes, the seasonal 
variation depending upon the orientation of the observing field. 
He accordingly believed the light-curve to be symmetrical, with a 
constant minimum of two hours’ duration due to the occurrence 
of an annular eclipse. That Yendell’s statement, however, may 
not be regarded as final is evidenced by the fact that Bemporad,° 
of the Osservatorio di Capodimonte, in 1914 published light-curves 
for both the June and December minima indicating a difference of 


t Astronomische Nachrichten, 109, 48, 1884. 2 Ibid., p. 59. 
3 Astronomical Journal, 9, 49, 1889. 

4 Astronomische Nachrichten, 157, 293, 1902. 

5 Astronomical Journal, 23, 213, 1903. 

® Astronomische Nachrichten, 199, 217, 1914. 
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o“25 in the depth of the curve. The June curve indicates also a 
secondary fluctuation during the minimum phase, and the Decem- 
ber curve suggests some similar peculiarity. 

Still more recently Shapley,’ in the discussion of a light-curve 
of U Cephei based upon Wendell’s extensive observations, recog- 
nizes a non-symmetrical tendency of the primary minimum and 
attributes the downward slant of the bottom of the light-curve 
to the irregular luminosity of the dark companion. He adds also 
that the close proximity of the components suggests a prolateness 
which might well show its effect upon the light-curve, although 
_ Wendell’s observations at maximum gave no definite indications 
of such a feature. 

With such discussions in mind, any new method of approach 
to the problem becomes of interest, especially if it affords oppor- 
tunity for measurements free from subjective errors. If any 
secondary fluctuation or other departures from the main curve do 
exist, it would seem that these could best be detected by a careful 
study of a single series of observations, since the very attempt to 
combine observations of different epochs and reducing to a mean 
curve would largely if not wholly mask any disturbance not com- 
mensurable with the principal period of the star. The advantage 
of photography in such a study is at once apparent, since it not 
only eliminates physiological difficulties, but makes possible a 
much larger number of independént observations during a given 
minimum than would be obtainable otherwise. With the thermo- 
electric photometer as a purely physical means for the measurement 
of the star-images, it would seem that any contributions through 
this means would be of special interest as a source of additional 
information of the behavior of the light-curve of this star. 

Plates taken for measurement by the method of diameters were 
found in the collection of the Yerkes Observatory, and a light-curve 
was published by Parkhurst in 1906.2 Plate R71, taken on 
June 25, 1905, with the 2-foot reflector, was used as the basis of 
the results then published. The same plate was the one first 
selected for measurement by the present method. Comparison 

* Contributions from the Princeton University Observatory, No. 3, p. 35, 1915. 

* Astrophysical Journal, 23, 79, 1906. 
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stars f (B.D. 81°30) and g (B.D. 81°27) were used. The exposures 
were made at approximately ten-minute intervals, the plate-carrier 
being moved by two revolutions of the screw after each exposure. 
Owing to the close proximity of the two series of images on the 
original plate, a second and somewhat enlarged negative was made 
from the original and used with the thermopile. 

The use of enlarged second negatives in this connection has 
the advantage not only of separating images in close proximity, 
but of increasing the contrast between star and plate ground, 
which may be reduced to clear glass if desired. This in itself 
greatly increases the precision with which thermopile readings may 
be made, since it affords much larger values of the galvanometer 
deflections and, with a fine-grained film, makes possible a more 
uniform background. On the other hand, it is doubtful if the 
enlarging process does not amplify errors in the original plate as 
well as introduce new ones in each reproduction. 


TABLE II 


PLATE R 71. JUNE 25, 1905. U CEPHEI 


Comparison Star B.D. 81°30, Magnitude 8.04 











Comparison Star B.D. 81°27, Magnitude 8.73 B=2.80 
l 
Image Mag. | Time Image Mag. Time 
| GMT. || G.M.T. 
ae 7.43 | ag*za%o@ |] t2......... 9.15 185 soto + 
See rege 7.88 43.4 a Seeorrrr 9.12 | 19 19.6 
ede dieidees 7.98 | 16 00.9 | T3----seeee 8.90 45.6 
Bee wacswe 8.07 20.5 P Ehecceseces 8.55 | 20 04.4 
acai hai eda 8.27 37.0 b BB wos esiecion 8.53 16.3 
_ oe 8.72 ge Se 8.23 28.2 
OS 9.14 | 17 20.8 f Seeerre 8.33 36.0 
GK arb in anes 9.30 39.7 | 18......... 8.21 44.0 
Doc ane ns 9.20 | 18 02.7 ig, OE ae 8.07 50.0 
eee 9.07 19.2 
t 











However, it is interesting to note that the resulting light-curve 
(Fig. 14) bears a striking similarity to that of Chandler to which 
reference has already been made. The times and magnitudes are 
given in Table II. The magnitudes of the comparison stars were 
taken as previously published, all points on the variable curve 
being referred to B.D. 81°30, of adopted magnitude 8.04. 




















THERMOPILE IN PHOTOGRAPHIC PHOTOMETRY 329 


Measurements were made on a number of other plates taken 
under widely differing circumstances, and indicated large residuals 
from the theoretical curve at minimum which could hardly be 
accounted for on grounds of errors in measurement. Repeated 
measuring gave the same results for a given plate with a surprising 
degree of accuracy. 

Through the courtesy of Professor E. C. Pickering of the 
Harvard College Observatory a number of plates, taken with the 

Mag. 


/ 





G.M.T. 15% 165 175 18 19" 20% 


Fic. 14.—U Cephei, minimum of June 25, 1905. Plate R 71 
x = points of half weight based on one comparison star 


“variable-star apparatus”’ attached to the 11-inch Draper telescope, 
were secured. By means of this attachment the variable was 
automatically photographed with a minute’s exposure at short 
intervals, the plate being displaced mechanically a fraction of a 
millimeter after each exposure. In this manner a long series of 
images may be impressed upon a 5 X8-inch (12}X 20cm) plate, 
covering a considerable portion of the minimum phase of the eclipse. 
The short exposure gave very small gray images of the variable 
when at minimum (9“2) which were difficult to measure. The 
photometer, however, was modified to meet the requirements, a 











33° HARLAN TRUE STETSON 


more efficient projection lens was substituted, a finding micro- 
scope was added above the stage, the galvanometer was removed 
to 2 m from the scale and read by a telescope magnifying 20 diam- 
eters. These improvements resulted in increasing the sensitiveness 
of the whole apparatus about five times and gave very good and 
consistent results. 

Plates C 11686 and C 11687 were taken during the minimum 
of December 29, 1898. From among the images available, those 
best suited for measurement were selected. Care was taken to 
compare the measurement of each star-image with the film back- 
ground in its own immediate vicinity, this precaution being ob- 
served to eliminate as far as possible any errors that might otherwise 
be introduced because of a changing transparency in the somewhat 
sky-fogged film. For the measurement of a single image not less 
than six settings were involved, the order of procedure being to 
set first on the image and then on the plate, alternating till the set 
of three pairs was secured. If the residuals appeared large for 
certain images, additional readings were sometimes taken in the 
order noted. The plate readings served also as a check on the 
voltage at the lamp terminals, and at the same time, when taken 
as indicated, eliminated any error introduced by the gradual 
decline, with consumption of current, of the potential of the storage 
battery. i 

As space would not permit the publication of the individual 
measures, a summary of the resulting data is given in Table III. 
The first column gives the number of the image on the plate, the 
second and third columns give the resulting A‘ for the variable 
and comparison star, respectively, and the fifth column contains 
the differences of the second and third columns multiplied by the 
plate constant 8, and therefore represents the differences of magni- 
tude between the variable and comparison star. The value of 
8 was determined from measures of stars B.D. 81°30 and B.D. 81°27, 
as was the case in the reduction of R 71. Columns 5 and 6 give 
the co-ordinates for plotting the light-curve, assuming the magni- 
tude of B.D. 81°30 as 8.04 as before. 

The graph resulting from these measures is shown in Fig. 15, and 
indicates a secondary fluctuation of about one-tenth of a magni- 
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tude, with an interval of 40 minutes between these apparent 
undulations. It is hard to account for this appearance on the 
grounds of systematic errors in measurement. The variable and 
comparison star were measured independently, and on different 
occasions, and it appears hardly likely that this curve is the deline- 
ation of accidental errors. Such a result might seem attributable 
to peculiarities in the plate or to atmospheric disturbances. Unfor- 
tunately, the overlapping of several rows of images on the plates 
restricted the use of suitable images other than those of comparison 
star B.D. 81°30. 


TABLE III 


Piates C 11686 AND C 11687. DATE OF PLATES, DECEMBER 29, 1898 


Comparison Star B.D. 81°30 





























at | 
Plate Image | Diff. 4 Mag. Mag G.M.T 
} | 
RS Kee 

C 11686 | I 0.750 | 1.054 0.304 I.220 | 9.260 10°59™23%9 
II 0.790 | 1.144 -354 1.436 | 9.476 | II 9 23.9 
21 0.830 | 1.191 .361 1.465 9.505 19 23.9 
31 0.799 | 1.122 - 332 1.345 9.385 29 23.9 
41 0.781 | 1.005 .324 1.315 9.355 39 23.9 
52 ©.790 1.136 . 346 1.404 | 9.444 5° 23.9 
62 0.801 1.162 . 361 1.465 9.505 I2 0 23.9 
72 0.845 1.174 | .329 1.338 9.378 IO 23.9 
83 | 0.837 1.163 | .326 1.322 9.362 2I 23.9 
Q2 0.815 I.150 | -335 1.360 | 9.400 3° 23.9 
102 | ©.740 a: . 382 1.551 9.591 4° 23.9 
112 | 0.781 1.140 | .359 | 1.457 | 9.407 5° 23.9 
122 | 0.760 a.118 | .358 1.453 9.493 13 00 23.9 
C 11687 2 0.833 1.116 . 283 1.145 9.188 9 23.8 
12 0.919 1.166 | . 247 1.003 | 9.043 19 23.8 
22 | 0.946 1.170 | .224 ©.909 8.949 29 23.8 
32 | 0.906 1.034 | .128 | 0.520] 8.560 39 23.8 
42 | 0.995 1.077 | .082] 0.333 | 8.373 49 23.8 
52 | 1.070 1.122 | .052 0.211 8.251 59 23.8 
62 | 1.056 1.068 | .O12 0.049 8.089 14 9 23.8 
fa 1.105 | 1.0907 | — .008| --0.032 | 8.008 19 23.8 
. 1.132 | 1.083 — .049 | —0.199 7.841 21 23.8 
| 93 | 1.144 | 1.090 | — .054 | —0.219 7.721 30 23.8 
1or | 1,164 1.118 | — .046 | —o.187 7.853 48 23.8 
112 1.125 I.105 | —0.020 | —0.081 7-959 59 23.8 





A set of three other plates was selected, Nos. C 9188, 9189, 9190, 
covering a portion of the decline, minimum, and rise of U Cephei 
under the date of June 13, 1896. The quality of the plates was 
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better than that of the previous pair, and it was possible to measure 
images at more frequent intervals. Two good comparison stars, 
B.D. 81°30 and 81°29, were used throughout. 

The results of the measures are given in Table IV. Columns 
1 and 2 give the number of the image on the plate, and the cor- 
responding Greenwich mean time of the exposure. Columns 3 and 
4 give the resulting magnitude of the variable as determined from 


Mag. 





G.M.T. 115 125 13> 14 15° 
Fic. 15.—U Cephei, minimum of December 29, 1898. Plates C 11686 and 
C 11687. 


the two comparison stars f and / respectively. The fifth column 
gives the direct difference of column 3 minus column 4 and forms 
the basis for suspected variability between f and h. The differ- 
ences of magnitude have all been referred to B.D. 81°30 as before. 
The magnitudes of B.D. 81°30 and B.D. 81°27 were adopted 
as for Plate R71, and the magnitude of B.D. 81°29 was deter- 
mined from the measures of the plates. 

The curve platted from these measures and represented in 
Fig. 16 was anything but what was expected, and indicated sources 
of difficulty not anticipated. Determinations of probable error 
pointed to variations of hundredths, whereas the curve advertised 
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discrepancies in the tenths’ place. Remeasures of offending points 
gave almost identical results with but slight variations in the last 
place, indicating that the source of trouble was not in the apparatus 


TABLE IV 


U CepHer, MIntmMuM FoR JUNE 13, 1896. J.D. 2413693 


. Comparison Stars Used: /, B.D. 81°30; g, B.D. 81°29; 4, B.D. 81°27 
Adopted Magnitudes: /, 8.04; g, 8.73; 4, 8.79 





SUMMARY PLATE 9188 
























































| Diff. Diff. 

Image G.M.T. | Mag.f. | Mag. h. | Mag. f- Image G.M.T. | Mag. f. | Mag. &. | Mag. f.- 
Mag. k Mag. & 

Gaede 14526™29%| 7.666 7.658 0.008 ee 15524™29*| 8.150 | 8.216 | —0.066 
BS vewet 28 29 7.704 7.678 .026 Pi sccs 26 29 | 8.112 8.048 064 
aay os | 32 29 7.670 7.764 .096 ices 29 29 | 3.170 8.332 |— .162 
Sse 34 20 7.718 7.768 |— .o50 ae 32 29 | 8.284 8.166 .118 
| ee 36 290 7.766 7.708 .058 ORis<<is 34 29 8.310 8.344 |— .034 
~ | 39 290 7.790 7.802 |— .o12 er 36 29 | 8.272 8.258 -O14 
ae 42 29 7.826 7.752 .O74 ee 39 290 8.286 8.418 |— .132 
ee 44 20 7.852 7.862 |— .o1o ae 42 29 | 8.430 8.488 |— .o58 
a 46 290 7.806 7.032 |— .036 io 44 29 | 8.468 8.530 |— .062 
See | 49 20 7.874 7.850 .024 aes 46 29 8.348 8.518 |— .170 
RO 52 29 7.862 7.804 |— .032 Sey 49 29 8.474 8.762 |— .288 
a. 54 29 7.800 7.824 |— .024 ae 52 29 | 8.392 8.356 .036 
OS 56 29 7.956 7.924 .032 ee 54 29 | 8.474 8.374 .100 
| AR 59 29 7.898 8.004 .094 "Se 56 29 | 8.678 8.680 |— .002 
Re man | 15 02 29 7.970 7.902 .068 ee 59 20 8.404 8.540 |— .046 
Oe uwe 04 29 8.138 7.802 .336 ee 1602 29 | 8.688 8.800 |— .112 
ee } 06 29 8.060 8.032 .028 || t00..... 04 29 | 8.560 8.634 |— .074 
eS | °9 29 8.130 | 7.988 -342 || 202..... 06 29 /| 8.788) 8.884 |— .0o96 
re I2 29 8.064 | 8.036 .oa8 || z0s..... | °o9 29 | 8.824 8.868 |— .044 
Oi ansal 14 29 8.086 7.97° a, ff. ae Ir 29 | 8.968 9.014 | .046 
"ee | 16 29 8.064 7.918 .146 | 109 ae 13 29 | 8.818 9.052 |—0.234 
ro Ties 19 29 | 8.064] 8.040] 0.024 | | 

en Se SEP ee } 

/ 

TABLE IV—Continued 

SUMMARY PLATE 9189 
| l l pry 
| Diff. Diff. 

Image | G.M.T. | Mag. f. | Mag. kh. | Mag. f- Image G.M.T. | Mag. f. | Mag. &. | Mag. f.- 

. | Mag. hk Mag. &. 

; | 

Bes ser | r6b32M29%/ 9.374 | 9.482 |—o. 108 2 Sr 17525™29*| 9.916 9.876 .040 

| Ria ie | 34 29 | 9.408 | 9.464 |— .056 . 28 29 9.658 9.662 |— .00o4 
Shiv s0x 37 290 | 9.206 9.416 |— .120 || 65..... 32 290 9.580 9.506 .074 

| | Se 39 20 9.282 | 9-464 |— .182 {| Beer 35 20 9.772 9.766 ‘ 

/ a i. ae L  Sceuanetinghanne 5 ae 38 29 9.764 9.768 |—0.004 
FEN Poe YN er See Fae | 7S+ee05 42 29 9.484 9.418 0.066 
iets. } 46 29 9.312 | 9.406 |— .094 || ae 45 20 9.538 9.448 .090 
as thes 50 29 9.388 9.400 |— .or2 |! 8z..... 48 29 9.618 9.508 .I10 
ies cen | 52 29 9.406 | 9.342 O65 1) BS... 006 52 29 9.488 9.348 .140 
ere 55 20 9.420 | 9.476 |\— .o56 || 88..... 55 20 9.466 9.348 .128 
OE  eeneies Ri ctlnh-ecoid Baaad steal } SP 58 29 9.710 9.530 . 180 
Sivesss 17 02 29 9.402 | 9 538 |— .136 re 18 02 29 9.630 9.428 . 202 

| ree 08 29 9.532 | 9.476 | O00 Tl Bosc. os 20 9.550 9.458 .092 
a 12 29 9.344 | 9.346 |— .002 || Too..... ek SVs Sota sta eie 6 Onwane 
a <chas I5 29 9.264 | 9.262 | .002 tea ok I2 29 9.460 9.388 072 
SBxivcws 18 29 9.320 | 9.418 |— .008 | 207 keves 14 29 9.554 9.438 °.126 
O6 Sccciak 22 29 | 9.588 | 9.748 |— .160 | | 
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TABLE IV—Concluded | 


SUMMARY PLATE 9190 





' = j i} l 








| Diff. || Diff 
Image | G.M.T. | Mag. f. | Mag. A. | Mag. f- || Image | G.M.T. | Mag. f. | Mag. h. | Mag. f.- 
| Mag. &. | Mag. &. 
Me as | 18h28™29"| 9.776 | 9.826 |—0.050 58.....| rgb22m2gs| 8.582 8.704 |—0.122 
Ee | 30 29 | 9.482) 9.426 |+ .056 , Pee 25 20 8.416 8.446 | —0.030 
ae hes | 32 29 | 9.506] 9.522 pes .016 || 63.....| 27 29 | 8.438 8.428 |+ .o10 
Rinses 35 20 | 9.338 | 9.440 in . 328 eee | 29 290 | 8.432 8.400 |+ .032 
eee 37 20) 9.388 9.444 |— .056 e's Wa | 32 29 8.538 8.536 |+ .002 
_. Veer 39 20 9-404 | 9.356 [+ Me Th 9B. as oe 35 20 8.350 8.370 |— .020 
Gi ass 42 20 9.254 | 9.380 |— .126 || 74..... OE Re A RPP ke 
Bases 45 290 9.2900 | 9.338 |— .048 || 77..... 4I 290 8.328 8.228 |+ .100 
— 47 20 | 9.528 | 9.556 |— .028 ae 45 20 8.302 8.142 |+ .160 
| ae 49 29 | 9.330 9.274 |+ .056 || 83..... 47 20 8.264 8.216 |+ .048 
~ SRR 52 29 | 9.222 | 9.322 |— .100 85.. 49 20 8.174 8.036 |+ .138 
SS ws 55 29 | 9.122 9.160 |— .038 aS 52 29 | 8.140 8.058 |+ .082 
Be ahs 57 29 | 9.182 | 9.318 |— .136 | 55 29 8.138 8.052 |+ .086 
y | RSS 59 29 8.778 | 8.868 |— .o90 || 93..... 57 20 8.196 8.084 |+ .112 
ee 19 02 29 | 8.918 | 9.010 |— .092 || 95..... 59 290 | 8.0904] 8.022 |+ .072 
is ones 05 29 | 8.742 | 8.770 |— .028 || o8..... | 2002 29/ 8.052 7.988 |+ .064 
Se > FF 9.070 | 9.136 |— .06 I0r.....| 5 29 8.090 8.038 |+ .052 
rae ONS Fae fl ] eae aR aan 07 29 7.984 8.040 |— .056 
Bes «3 I2 29 8.712 | 8.856 |— .144 !| 105..... | °9 20 8.038 7.828 |+ .210 
PS 15 29 8.608 | 8. — .082 |} 108.....| I2 29 7.932 7.884 |+ .048 
awed 17 29 8.764 8.724 |+ .040 || I10..... | 14 29 7.920 7.824 |+0.096 
EES. | 19 29 | 8.876 | 8.748 |+ 








nor in the method of measurement. The discrepancies were there- | 
fore photographic, to be interpreted as irregularities in the exposure | 
due to-rapid changes of atmospheric transparency, as an uneven | 
condition of sensitiveness in the plate film, or, as seemed hardly 
likely, due to actual changes of light among the stars in question. 

The singular feature is that for the most part the general char- 
acter of the light-curve is the same for either comparison star, 
indicating that the chief cause of the discrepancy is in the images | 
of U Cephei. If the large irregularities are due to atmospheric 
changes, it is difficult to see why the variable alone should be 
affected, while the two comparison stars fail to show a corresponding 
disturbance. The most outstanding difficulty is the behavior of 
the curve during minimum, around 1720", where, within the short 
space of half an hour, both comparison stars indicate the variable 
to experience a change of nearly half a magnitude in light. It has | 
seemed wiser in this preliminary graph to leave the points connected 
with straight lines rather than.to bias judgment by any attempt at 
drawing a mean curve until more complete data are at hand. The 
combination of several such curves at minimum would doubtless 
“smooth out”? many irregularities, but the argument perhaps 
needs renewed emphasis that the combination of phenomena at 
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separate minima only tends to mask any singularities not a definite 
function of the 2.49 days. 

Another fact at once apparent from Fig. 16 is that the curves 
drawn from the data furnished by the two comparison stars cross 
at stated intervals. The curve belonging to B.D. 81°30, for 
example, clearly lies above that drawn for B.D. 81°29, from 165 
to 17, whereas from 17530™ to 18%30™ the reverse is true. A 





G.M.T. _15® 165 17% 184 19* 208 


Fic. 16.—U Cephei, minimum of June 13, 1896 
@ = comparison star f; © = comparison star h 


third crossing of the two curves occurs at about 19°30". This 
behavior suggests a variability in one or the other of the comparison 
stars, and becomes more apparent when we examine the two com- 
parison stars platted against each other as in Fig. 17. Here, in 
spite of the large residuals from a mean curve, the rise and fall of 
the curve at an approximate interval of an hour and three-quarters 
is at once evident. Sufficient measures of a third star, B.D. 81°26, 
were made to indicate that the eighth-magnitude star B.D. 81°30 
was the cause of the discrepancy. We must therefore either try 
to account for the rise and fall in Fig. 17 on the ground of a progres- 
sive undulating change in atmospheric transparency continuing 








Diff. Mag. f —Mag. h 
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throughout the entire exposure, but not appreciably affecting the 
adjacent star, or else suspect the comparison star B.D. 81°30 of 
variability. If this should be the case, it might well be a contribut- 
ing factor in the cause of the various forms of the light-curve of 
U Cephei previously published, when this star has been the principal 
comparison star used." 

Whatever may be gleaned from these preliminary results, it is 
hard to see how a straight flat-bottom curve, as demanded on the 
grounds of an eclipse theory alone, can be reconciled with the obser- 
vations. Ifthe binary system be held to account for the main period 


+0.5 : 





—0.5 
G.M.T. 15% 165 175 18 19 20h 21 


Fic. 17—Curve showing suspected variation in B.D. 81°30 


of the star,? it may seem necessary to suppose other vital factors 
operative in the light-curve of the star. In any event, enough 
has been exhibited already to indicate that with every determina- 
tion complexities arise which forbid an easy verification of anything 
more than a mean light-curve of general characteristics. With the 
continuance of the method of investigation here outlined it is 
hoped that the completion of further measures in all parts of the 
light-curve may at length afford sufficient data for the formation 
of more conclusive ideas regarding the behavior of the light-curve 
of this most interesting object. 

* Compare with the results of Capodimonte, to which reference has already been 
made (Astronomische Nachrichten, 199, 217, 1914). 


2 For a discussion of irregularities in the main period, see Yendell’s article (Joc. cit.). 
Further investigation of this point is being undertaken by Shapley. 
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In order to determine the order of precision which could be 
expected in the measurements of the Harvard plates, the probable 
error was computed from a number of measures of the same star- 
image, of U Cephei near minimum, and found to be +o™“oog for 
twenty: settings or corresponding to about +o™“o25 for a set of 
readings as ordinarily taken for each star. Results from measuring 
diameters with a micrometer microscope and reducing to magni- 
tudes by means of the square-root formula yielded a probable 
error nearly five times as great. This large value resulted from 
the fact that, the plate being unguided, the images were somewhat 
enlarged, and the variable was so faint at the minimum that the 
images then taken were still gray and showed no well-defined 
periphery. The use of the thermopile promises well for the 
measurement of star-images not suitable for reduction by the 
diameter method. With longer exposures on the variable at mini- 
mum and special care in the development of the plate it should 
be possible to reduce the probable error to a still lower figure. 
Experiments with the ferrous-oxalate developer show a great ad- 
vantage to be gained by its use for photometric purposes, especially 
for plates to be measured thermo-electrically, inasmuch as good 
contrast results without the slightest possibility of chemically 
fogging the unexposed portion of the film. 

Sources of error in the film itself are without doubt a cause of 
many discrepancies, and further refinement must be looked for 
in adopting a modus operandi which will reduce the effect of such 
residual errors to a minimum. 


CONCLUSION 


In the development of the present method of investigation, the 
first consideration has been to adapt a piece of apparatus to the 
measurement of stellar magnitudes from photographic plates which 
should eliminate so far as possible some of the difficulties incident 
to other methods. In so doing care has been taken to consider the 
errors involved and to reduce them to a minimum. The degree 
of precision attainable by means of the thermo-electric measure- 
ments has been tested under widely different circumstances and 
enough numerical results tabulated to indicate what additional 
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contributions are to be expected from the use of such an instrument 
in stellar photometry. 

The several lines of special investigation suggested by the results 
of these preliminary determinations indicate a fruitful field for 
extended study in many directions. It has seemed wise, therefore, 
in the present treatment to outline in considerable detail the 
essential steps in the development of the instrument and methods 
of reduction, with a preliminary application to several different 
problems, rather than to devote an early investigation in a new field 
to a more exhaustive treatment of a single topic. 

The precise determination of magnitudes by means of the 
thermo-electric method of measurement has already been discussed. 
The difficulties mentioned are not new to the present method; 
rather, they are those which have long occupied the attention of 
investigators in stellar photometry. It is to be hoped that the 
measurement of plates in the manner herein described, by reducing 
another element of uncertainty, may facilitate more rapid progress 
in the field of photometry. 

As has been shown, the thermopile may be used for the direct 
measurement of plate opacities, covering the same field as that of 
the Hartmann microphotometer with about the same degree of 
accuracy. The points of difference between the two instruments 
when so used may again be mentioned. The reductions of the 
Hartmann readings depend upon the calibration of the photographic 
wedge used in comparison, whereas with the thermopile arrange- 
ment the deflections of the galvanometer may be taken at once as 
proportional to opacities on the plate. The thermopile further 
has the advantages of restricting a much smaller area for examina- 
tion than the microphotometer, and of making the comparison 
entirely independent of physiological effects or personal equation. 

In measuring opacities of extended areas such as are found on 
solar, lunar, or nebular photographs, these advantages may be 
in part offset by the uncertainty in the assumption of constant 
absorption of glass and gelatine throughout the plate. Such a 
source of error will not prove troublesome whenever adjacent 
settings on the comparatively uniform unexposed portion of the 
film can be made. 
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The ease with which the apparatus may be adapted to the service 
of spectrophotometry in the measurement of objective-prism plates 
affords opportunity for an important study of the relation of color- 
index to spectral type, and may afford further means for a more 
precise definition of ‘‘magnitude.”’ 

In a similar manner the application of such spectral measures 
to the determination of the radial velocities of faint stars photo- 
graphed with the objective prism is a promising field for investiga- 
tion. Insucha case it should be possible to locate with considerable 
accuracy the position of the absorption band of neodymium used 
as a comparison spectrum. The method becomes equally appli- 
cable to line-of-sight measures for stars having broad-line spectra. 

A line of work most promising for the thermo-electric photom- 
eter is that of the photographic study of short-period variables. 
In the case of fields where close comparison stars are available 
and a long series of exposures are made on the same plate, the large 
and uncertain errors of atmospheric transparency, plate con- 
stants, and reduction to a standard scale are reduced to a mini- 
mum and the accuracy of the measures appears to be limited only 
by the quality and uniformity of the plate. Combining the ad- 
vantages of photography with refinement of measurement, the 
thermopile so used may yet be able to accomplish for the fainter 
stars what the photo-electric cell and the thermopile in a more 
direct way are promising for the stars of brighter magnitudes." 

In conclusion, the writer wishes to express his appreciation to 
both astronomers and physicists who in any way have aided in 
the progress of this investigation. He especially wishes to acknowl- 
edge his indebtedness to Professor Fox, of the Dearborn Observ- 
atory, Northwestern University, for the provision and loan 
of important pieces of apparatus; to Professors Michelson and 
Millikan, of the University of Chicago, for so kindly placing at 
his disposal the resources of the Ryerson Physical Laboratory; to, 
Professor E. C. Pickering, for the loan of valuable plates from the 
Harvard College Observatory, and for many important suggestions; 
to Professor Frost, for the use of the facilities of the Yerkes 


1 See W. W. Coblentz, “Radiometric Measurements of 110 Stars with the Crossley 
Reflector,” Lick Observatory Bulletin, 8, 104, 1914. 
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Observatory, and for his kindly interest in the progress of the work; 
and to Professor Parkhurst, for the use of certain plates taken by 
him, and for that constant appreciation and material encouragement 
which have done much to further the application of laboratory 
methods to the problems of astronomical photometry. 


SUMMARY 


The results of the foregoing paper may be briefly summarized 
as follows: 

1. An instrument has been devised, using a delicate thermopile 
as a receiver, whereby the magnitudes of stars may be determined 
from either focal or extrafocal plates by purely physical means. 

2. A reduction formula has been adopted which shows a fourth- 
root relation between galvanometer deflections and stellar magni- 
tudes, and this relation has been checked by the measurements of 
magnitudes from plates taken under various conditions, and with 
four different telescopes ranging from six to forty inches in aperture. 

3. Provision has been made for the adaptation of the apparatus 
to the measurement of spectral intensities, and the method illus- 
trated in the measurement of objective-prism spectra of B-type 
stars for purposes of spectral photometry. 

4. An extensive study of the eclipsing variable U Cephei has 
been begun, the preliminary measures here presented indicating 
the presence of light-changes not explained on the eclipse theory 
alone, and revealing a suspected variability in the much-used com- 
parison star B.D. 81°30. 

5. From the measurements of the variable-star plates a some- 
what higher degree of accuracy is to be expected from thermo- 
electric measures than from measures of diameters, the thermopile 
promising distinct advantage in the case of images too ill-defined 
to be suitable for the micrometer microscope. The principal 
sources of error appear to be the changing atmospheric conditions 
during the exposure at the telescope and irregularities in the 
photographic plate. ; 


YERKES OBSERVATORY 
November 1915 











THE PRODUCTION IN THE ELECTRIC FURNACE OF 
THE BANDED SPECTRA ASCRIBED TO TITANIUM 
OXIDE, MAGNESIUM HYDRIDE, AND CALCIUM 
HYDRIDE' 

By ARTHUR S. KING 

The purpose of this investigation was to produce the banded 
spectra of titanium oxide, magnesium hydride, and calcium hydride 
in the tube-resistance furnace, to compare the conditions required 
for their appearance with those observed in the arc and spark, and to 
note any phenomena resulting from the special control of conditions 
possible with the furnace. 

All of these banded spectra have been identified in the spectra of 
sunspots, and the bands attributed to titanium oxide also occur 
in the spectra of third-type, or Antarian, stars, which fact has led 
Fowler to designate them as the “Antarian bands.” ‘The presence 
of such banded spectra is taken as evidence of a relatively low 
temperature, while the appearance of the titanium oxide bands 
furnishes direct evidence of the presence of oxygen in third-type 
stars and in the sun. 


THE TITANIUM OXIDE BANDS 


An extended investigation by Fowler? furnished concordant 
evidence that the Antarian bands in the arc and spark result 
from the presence of titanium oxide. The flutings occur in a series 
of groups, all degraded toward the red, throughout most of the 
visible spectrum. Fowler gives the following wave-lengths for the 
first head in each group: 


4353.68 5597.92 
4462.34 5760.15 
4584.62 5954.66 
4761.08 6158.86 
4954.78 6357.9 
5167.00 6651.5 
5448.48 7054-5 


* Contributions from the Mt. Wilson Solar Observatory, No. 114. 
* Proceedings of the Royal Society, A, 79, 509, 1907. 
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It was noted by the writer,’ during an examination of the line spec- 
trum of titanium in the furnace at different temperatures, that these 
bands did not appear with the furnace im vacuo. Experiments 
were recently undertaken to see whether the use of oxygen would 
bring them out. The photographs were made with a 15-foot 
concave grating spectrograph, this being set for the first observa- 
tions to cover the region from \ 5500 to X 7200 in the first order, 
with a dispersion of 1mm=3.7A. The furnace was usually 
operated at about 2300° C., which gives a strong titanium spectrum 
without excessive vaporization of carbon. 

Admitting air to one-half atmospheric pressure gave no bands. 
They again failed to appear when the furnace was full of air, the 
outlet valve being left open. It seemed probable that a strong 
current of oxygen was needed, since, there being no circulation 
through the horizontal tube, the incandescent carbon probably 
used up the small supply of oxygen before the titanium was fully 
vaporized. The windows in the ends of the furnace chamber 
were therefore removed, a silica tube inserted in one end of the 
furnace tube, and a strong blast of air forced through by means of 
a bellows. The whole series of bands then came out with great 
brilliancy. Temperatures of 2300° and 2600° C. were used with 
no perceptible difference except in exposure time. A noteworthy 
effect was the almost complete suppression of the line spectrum of 
titanium when a strong current of air was used, though the lines 
of several impurities appeared as usual. Evidently a sufficient 
supply of oxygen resulted in turning all of the titanium vapor into 
the oxide, which seems to give only the band spectrum. 

Pure oxygen in the furnace was next used, this being led directly 
into the tube from a tank of compressed gas. The spectrum could 
be changed at will from a mixture of lines and bands to a pure 
band spectrum by increasing the current of oxygen. When the 
bands were at their strongest, the image on the slit was a greenish 
yellow, similar to the envelope of a long carbon arc containing 
titanium, which is one of the best of the usual sources for obtaining 
the band spectrum. 

* Mt. Wilson Conir., No. 76; Astrophysical Journal, 39, 139, 1914. 
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There was no evidence of a material change in temperature 
caused by the introduction of oxygen. The same exposure time 
with and without. oxygen gave the same general intensity of 
spectrum, this being composed of bands in the one case and 
lines in the other. Pyrometer measurements of the interior of the 
tube were practically the same in the two cases. Allowing for 
differences in clearness of the vapor, the change of temperature 
could scarcely have been greater than 50°. 

Comparing this source with the usual flames employed in produ 
cing spectra, we have here the metallic vapor heated by the graphite 
tube instead of by a union of flame gases, and then the process of 
oxidation taking place. The results are probably essentially the 
same as in a flame of the same temperature, but the furnace with 
oxygen gives a method of obtaining a series of flame temperatures 
for which different combinations of gases are usually required. 
There is also the possibility, by the use of very high furnace tempera- 
tures combined with a stream of oxygen, of exceeding any of the 
flame temperatures regularly attainable. 

It may be noted that a distinct difference in gradation of 
intensity with the wave-length appeared for the titanium band 
spectrum in the furnace as compared with the arc, the bands at 
the red end being relatively stronger in the furnace. 

A supplementary set of photographs was made for the region as 
far as \ 4200 with the furnace in vacuum, in air at atmospheric 
pressure, and with varying amounts of oxygen. The results ob- 
tained for the yellow and red region were repeated. By the use of 
a strong current of oxygen, so that a flame extended several centi- 
meters from the end of the furnace tube, a strong band spectrum 
appeared, while the titanium line spectrum was eliminated except 
for the faint appearance of a few of the lines strongest in the 
regular furnace. It was noteworthy, however, that the more promi- 
nent lines of iron and chromium, as well as \ 4227 of calcium and 
\ 4607 of strontium, appeared as strongly as could be expected 
from the slight impurities present. The evidence thus seems fairly 
conclusive that the bands of titanium oxide have been correctly 
ascribed to an actual compound, since otherwise, without a large 
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reduction of temperature, there would seem to be no reason for 
the disappearance of the line spectrum. 

Reproductions are given in Plate VI of the titanium spectrum 
from \ 4427 to 5065 in which the four strips represent respectively 
the arc in air; the furnace with a moderate stream of oxygen, both 
lines and bands appearing; the furnace with a strong stream of 
oxygen, the spectrum consisting almost wholly of bands; and the 
furnace in air without extra supply of oxygen, in which case the 
bands are absent. The wave-lengths are those given by Fowler 
for the principal heads. 


THE BANDS OF MAGNESIUM IN HYDROGEN 


In a paper giving detailed measurements of the fluted spectrum 
obtained when magnesium is burned in an atmosphere of hydrogen, 
Fowler’ summarizes previous investigations as to the dependence of 
the spectrum on hydrogen and as to the probability of its corre- 
sponding to a chemical compound. The bands consist of three 
main groups, beginning at AX 4844.9, 5211.1, 5621.6 and shading 
toward the violet. They appear best with the arc or spark in an 
atmosphere of hydrogen, but may often be seen when only a small 
residue of hydrogen is present. The furnace results are in close 
agreement with those for the arc, which is a fact of some signifi- 
cance, since the arc discharge is greatly modified by a hydrogen 
atmosphere, and this might account for the appearance of another 
spectrum. 

The strongest of the bands, with its head at \ 5211, was situated 
most favorably for observation in the furnace photographs, having 
a portion of its length undisturbed by the carbon spectrum. The 
first experiments were made by simply admitting hydrogen to the 
furnace chamber. Photographs at 1900°, 2200°, and 2500° C. were 
made for each of three pressure conditions, viz., with a pressure 
below 1 cm of mercury, the residual gas being air, with hydrogen 
at a pressure of 10 cm, and with hydrogen at atmospheric pressure. 
The results were somewhat variable, the main effect being that a 
high temperature was unfavorable for the bands, which, if they 
appeared at 2500, were usually in absorption. The band at 


* Philosophical Transactions of the Royal Society. A, 209, 447, 1909. 
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dX 5211 occasionally appeared distinctly with the furnace pumped 
out. The experiments with the tube-arc have shown, however, 
that considerable hydrogen is occluded by the graphite tube and 
probably by the metallic magnesium, since hydrogen linés appear 
strongly in the tube-arc under these conditions. 

It appeared possible that the effect of hydrogen had not been 
fairly tested, since a certain degree of combination with the carbon 
would occur, probably with the formation of acetylene. Arrange- 
ments were therefore made for leading a stream of hydrogen directly 
into the furnace tube by means of a silica tube passing through a 
rubber stopper in one end of the furnace chamber. The hydrogen 
was generated chemically from zinc and sulphuric acid and dried 
by concentrated sulphuric acid and sodium hydroxide. Passing the 
gas through the tube in this way, the pressure being held at 1o- 
20 cm, proved very effectual in producing the bands. They were 
obtained in emission at 2500°, but much better at 2200°, when the 
band at A 5211 was stronger than the carbon head at \ 5165. The 
line spectrum of magnesium appeared with the bands, though the 
b group and \ 4571 are the only lines given strongly by the furnace 
in this region. 

The experiments seemed to indicate a clear dependence of the 
bands on the presence of hydrogen, without any apparent change 
in the action of the source. 


THE BANDS OF CALCIUM IN HYDROGEN 


The calcium arc burning in hydrogen was observed by Olmsted" 
to show a banded spectrum in the red with well-marked heads at 
dA 6382.2 and 6389.3 and a fainter group in the B region, in both 
cases shaded toward the violet. The furnace was supplied with a 
stream of hydrogen through the tube, as in the experiments with 
magnesium, and the calcium bands were found to be very similar in 
behavior to those of magnesium. They frequently appear, as 
was noted by the writer,’ with the furnace pumped out, apparently 
from occluded hydrogen. A current of hydrogen through the tube 
brings them out strongly at temperatures from 1900° to 2200”. 

* Mt. Wilson Contr., No. 21; Astrophysical Journal, 27, 66, 1908. 

2 Mt. Wilson Conir., No. 35; Astrophysical Journal, 29, 190, 1909. 
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Higher temperatures are not so favorable, though the bands have 
been obtained at 2500°. 

The component lines of the band to the red of \ 6700 can be 
measured to advantage when furnace photographs are made with 
higher dispersion. The denser portion shows distinct maxima at 
approximately \ 6902 and \ 6919, while a more scattered structure 
extends farther to the red. The calcium line spectrum appeared 
regularly with the bands. 


SUMMARY 


1. The bands ascribed to titanium oxide show a clear depend- 
ence on the presence of oxygen in the furnace, the spectrum con- 
sisting of bands alone when a sufficient amount of oxygen is supplied. 

2. The bands attributed to magnesium hydride and to calcium 
hydride require in each case the presence of hydrogen in the furnace 
to give the bands strongly, although they frequently appear with 
moderate strength when only a small quantity of hydrogen is 
present. 

3. While the bands appear through a considerable range of 
furnace temperatures, about 2300° C. seems to be the upper limit 
for their greatest strength. mart 


Mr. WItson SOLAR OBSERVATORY 
April 5, 1916 























RESEARCHES ON SOLAR VORTICES' 
By CARL STORMER? 


I, INTRODUCTION 


Professor Hale, as is well known, has done fundamental work 
on sun-spots in discovering the Zeeman effect in their spectra, thus 
proving the existence of very strong magnetic fields in the spots 
and in the surrounding whirl revealed by the spectroheliograms. 
In an attempt to employ mathematical methods in the investi- 
gation of these phenomena, I spent some months of the summer of 
1912 at the Mount Wilson Solar Observatory. 

The detailed study of the collection of spectroheliograms made 
at the Observatory suggested to me an application of the classical 
researches on terrestrial cyclones made in 1876 by my countrymen, 
Guldberg and Mohn.’ The hydrogen flocculi seemed to be arranged 
around a single spot in curved paths that were very similar to 
logarithmic spirals; and just the same curves were found by Guld- 
berg and Mohn as trajectories of the air particles in the outer part 
of cyclones. It therefore seemed advisable to start with the 
hypothesis that the motion of charged electric gas molecules takes 
place along such spirals around the sun-spot center, and compute 
the resulting magnetic field. 

This investigation which I made in Pasadena showed that the 
lines of magnetic force due to a plane whirl of the kind mentioned 
above are space-curves whose projection on the plane of the whirl 
are also logarithmic spirals intersecting the first at right angles. 

This result led me to the idea first advanced by Brester‘ in 
1909, that the hydrogen flocculi resemble terrestrial auroras, thus 


* Contributions from the Mount Wilson Solar Observatory, No. 109. 

2 Research associate of the Carnegie Institution of Washington, Mount Wilson 
Solar Observatory. 

3“ Etudes sur les mouvements de V’atmosphére, Premiére partie’; Programme 
de l'Université pour le 2° semestre, 1876, Christiania, 1876; and Meteorologische Zeit- 
schrift, XII. Band, 1877, No. 14. 

4 “The Solar Vortices of Hale,’ Proceedings of the Amsterdam Academy of Sciences, 
January 30, 1909. 
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implying that the visible whirls are not real current-lines, but lines 
of magnetic force due to a magnetic field at a lower level—a view 
expressed by Deslandres‘ in 1910. Professor Hale had previously 
stated in his papers? that the Zeeman effect was probably due to 
a low-level vortex. In accordance with the mathematical result 
mentioned above, the current-lines of this hypothetical whirl at 
a low level would be logarithmic spirals tending more and more to 
become circles according as the structure of the hydrogen whirl 
tended to be radial. This would much better explain the strength 
of the Zeeman effect of sun-spots of radial structure. 

In the following pages is given a detailed account of my 
researches on the mathematical theory of solar vortices. 


II. DEFINITION OF THE WHIRL 


Let us first consider a logarithmic spiral in the XY-plane of a 
system of rectangular co-ordinates in space. If we use polar 
co-ordinates, R and ¢, in that plane (see Fig. 1), the equation of 
the spiral is 

R=Ae? tan w 
where A and w are constants. A only fixes the position of the 
spiral; for if we introduce 


$=—cotwln A 


we get 
A=e —®, tan w 


and consequently 
R=e(—%.) tan o 


which is the spiral 
R= e? tan °. 


turned through the angle ¢, around the Z-axis. 

The other constant, w, determines the shape of the spiral. In 
order to get all possible values of tan w between — « and + ~, we 
can suppose w to be within the interval —z/2 to +7/2. 

* Comptes Rendus, January 10,1910. See also G. Hale, “Notes on Solar Magnetic 
Fields and Related Phenomena,” Publications of the Astronomical Society of the 
Pacific, 22, 71, 1910. 

2 Annual Report of the Solar Observatory of the Carnegie Institution of Washington, 
1909, p. 166 (24). 
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On the other hand, w is the constant angle between the tangent 
MT and the normal MN on the radius vector, both reckoned in the 
direction of the increasing angle ¢. w is negative if MT lies 
between MN and the origin, positive if MT lies on the other side. 
If w is negative, R will decrease with increasing ¢, and we approach 
nearer and nearer to the origin in following the spiral. If w is 
positive, R will increase with increasing ¢, so that we recede from 
the origin in following the spiral in the direction of increasing ¢. 

In Fig. 3 are shown various types of spirals corresponding to 
various values of w. Each separate figure contains a series of 
congruent spirals seen from the positive side of the Z-axis, with 
the origin in the center; the corresponding value of w is written 


Z 











? N 
Y * . 


Fic. 1 Fic. 2 


above the figure. w=o corresponds to circles concentric around 
the origin. If w approaches —x/2 or +7/2 the spirals tend to 
become straight lines through the origin. 

After this review of the well-known properties of logarithmic 
spirals, we will consider an ideal vortex defined as follows: 

The whirl shall be bounded by two cylinders of revolution 
with the Z-axis as axis and with radii equal to p, and p., by the 
XY-plane, and by a parallel plane at the distance # above it. We 
consider h small as compared with p, in order to have a flat whirl. 
The motion of the gas in the space mentioned is supposed to be 
permanent, and the trajectories of the gas particles are supposed 
to be logarithmic spirals parallel to the X Y-plane, whose projections 


on that plane are 
R= Ae? tan w 
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where w is the same for all particles, and where A is determined 
by the position of the-spiral. Let us suppose further that the 
velocity and electric charge in unit volume are functions of the dis- 
tance from the Z-axis only. 
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It is clear that such a whirl does not exist in reality; there must 
be a source and a sink for the coming and going masses of gas in the 
whirl, but in the first approximation we will neglect the electric 
action of masses outside this hypothetical whirl; their action will 
be calculated later, when we have discussed the formulae for this 
first simple case. 
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To fix the motion without ambiguity, we will call the velocity 
v positive if the motion corresponds to increasing ¢, negative if it 
corresponds to decreasing ¢. Thus the motion seen from the 
positive side of the Z-axis will be clockwise for positive v, counter- 
clockwise for negative v; and the gas will approach the Z-axis if 
the product wv is negative, and recede from the Z-axis if it is 
positive, as is seen in Fig. 4. 
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Let us now calculate the magnetic action of this hypothetical 
whirl on a unit of north magnetism placed at a point outside the 
whirl. For the sake of brevity we first suppose w, v, and the elec- 
tric charge to be positive. In order to have a solenoid element 
bounded by current-lines, let us consider an element of the whirl, 
bounded by the following surfaces: 

The two planes 
s={, 2={+¢; 


the two cylinders 
R=p, R=p+Ap; 
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and the two cylinders 


R=Ae? tan o. R= Ae (et) tan °, 


where Af, Ap, and Aé@ are positive infinitesimal quantities, and 
Aé and A¢ infinitesimal compared with Ap. 
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Let us make use of the c.g.s. system of units, and the electro- 
static system for the electric charge. The cross-section of our 
element will be 

pA sin w A. 


If € is the electric charge in a unit volume of the gas, there will 
pass as a convection current through this section in one second 


evp sin w AOAé electrostatic units, 
and consequently 


Pan _ M . . 

3X10" evp sin w AGA electromagnetic units. 
If by i we denote the current in amperes passing through the 
element, we have 


I 
i= ——— evp sin w AOAL. 
3X10" . ‘ 
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On the other hand, the length of the element is 
As= oe ; 
sin 


Now, as is well known, the magnetic action of this current-element 
on one unit of north magnetism situated at the point M with 
co-ordinates x, y, and z will have the following components parallel 
with X-, Y-, and Z-axes: 


1 .As 
p:=— i — (cos B, cos y.—cos y; cos B,) , 

Io 9 

1 .As 
py=— i; (cos y; COS a,—COS a; COS y2) , 

I 


A 
p= —i ; (cos a, cos B,—cos B, cos a) , 
where 7, is the distance from the element to the point M; az, 8:, ¥:, 
the angles between the current-element and the positive directions 
of the X-, Y-, and Z-axes; and az, B., 72, the angles between the 
direction from (M) to the current-element and the same three axes. 
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Now choose the point M in the XZ-plane (Fig. 6) so that: 
y=o. 
The co-ordinates of the current-element are 


pcos@, psin@, f, 
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and consequently 


=(x—p cos 8)?+(p sin 6)?+(2—¢)? 
and 
p cos 0—x p sin 6 c—% 


a. ——— cos = cos eam 
2 T; c B; rT; , * rT; 





On the other hand, we have 
a,=—— 06, 8, =6—o, we 4 
which give 
cos a,=sin (w—@), cos B,=cos (w—8), cos y:=o0. 


If we substitute the values of the angles a and 8, we get 


pea 2 iO 005 (w 6) 
A , 
py= SiS ERO cin (w—6) 
= +4 3 {x cos (w—8)—p cos w}. 


These formulae are valid for ©, v, and w positive. But if 
we change the sign of © or of v, the magnetic action will change 
direction, that is, p,, ~,, and p, will change their signs. But this 
will also be the case with 7 and consequently the formulae are 
valid for both signs of the quantities v and ¢«. In an analogous 
manner we see that they are valid for all values of the angle w. 

The components of action of an element situated symmetrically 
with respect to the XZ-plane are obtained by changing 6 into —6; 
we find for these components 





»,__ 1 ,As (s—Q) 3 
pr= ae r cos ( +6), 
“ae jas oP 2—€) i 
 - P sin (w+) , 


. As 
p= is {x cos (w+6)—p cos w}. 


I 
Io 
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If we introduce the values of 7 and of As and add, we get 





























,_ 2COS® ApAlA@ 
pit+p.= 3x10" ” (s—{) epv cos 6 ma aa 
,_ 2S8ino ApAlAé 
Pyt p= 3X0" * (s—{) epv cos 0 s 
,_ 20S ApALlAé 
Pst+p.= 3X10" (x cos 6—p) ep r 
From these the action of the part of the whirl bounded by the 
two planes sf, s={+Al, 
and the two circular cylinders 
R=p, R=p+Ap 
is obtained by an integration with respect to 6, which gives the 
components (rm 
2 COS w cos 6 dé 
ain idee ApA — 
3X10" io 0) een See > 
eo 
2 sin e cos 6 dé 
i) 
iB... Rael ON ApA ec web. 
3xro~ @— 6) £90 Apt no? 
e/o 
2 COS w " & Cos 6—p 
ApA 6. 
3X10” epv Op at Pp d 


Here e and 2 are independent of the angle @ according to the suppo- 
sition mentioned above that they are functions of p alone. 

The action of a ring bounded by the same two planes and the 
two cylinders R=p, and R=p, is obtained by an integration with 
respect to p, which gives the components 


2 COS w 


3x ro’ 





2 sin w 
3x 10” 


2 COS w 
3X 10” 





(* Po (?™ 
(z—{) Ag n( 
e/py eo 
(* Pz (* 9 
(z—{) A¢ cn 
e/ 


eV) 








Pr 


Here e and » are functions of p only. 








Po ite 
a | £pv (f stot ae) 
Py ° ; 
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A final integration with respect to ¢ gives the components of 
the action of the whirl itself: 


h p T 
2 COs w . cos 6 dé 
~ 3X10" (2—¢) | Epv ({ te) | dg 
° Py ° 


with analogous expressions for the other components. 

But here we suppose the product epv to be a continuous function 
of p; and the point (x,0,z) being outside the whirl, 7, is positive 
over the entire field of integration. This field being also finite, 
the succession of the integrations is indifferent, so that we can 
write more simply 


2 COS w ; a . cos 6 
vie -— ” dOdpdt. 
 eicuill { { { ; ° ed " ; _ 
° I ° 


On the other hand, as everything is symmetrical around the 
Z-axis, we get the components in an arbitrary point (x,y,z) by 
substituting R= V x*+-y* for x; hence 








Ph (oe, (2 
(z—£) epv cos 6 



































2 COS w 
a0 —-———— — d0dpdé 
3X10 soe Sap Ds 
, Ch Pe, (em 
2 sin w (s—{) epv cos 0 
Hs= d6dpdé (1) 
10 3 
3X10 eJ/o e/py eo D 
2 COS w Dee! Sg Neg (R cos 6—p) epv 
H, a * 10” is ae. d6dpdé 
3 eo ep, e/o 
where 


D?= R?—2Rp cos 0+-p?+ (s—{)’. 


Hr is the component normal to the Z-axis, H, the component 
parallel to that axis, and H, the component normal to Hz and H3. 
As regards the signs, Hp is positive in the direction from the 
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Z-axis, H, is positive in the direction of increasing z, and Hy posi- 
_ tive in the direction of increasing ¢ (see Fig. 7). 


Z ,H, 











4 
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III. PROPERTIES OF THE COMPONENTS OF THE MAGNETIC ACTION 
EQUATIONS OF THE LINES OF MAGNETIC FORCE IN SPACE 


Let us, in the formula for Hr, Hy, and H,, introduce the function 


Epv cos O sae , 
tof Cfertine ot 


As the point (x,y,z) is outside the field of integration, and 





oe 5® 
consequently D>o, we obtain the derivatives = and % 35 follows: 




















bR 
(*h (*p 4 (*™ 7 
8 2 ° 8 R cos 8 
8R 3X10" <P? | 8R 5 
ev e/p; 4 eo ‘ 
sb 2 & gar. 3 "" Rcos 8 a6) 1 ae 
bs ~ 3X10" pn | 8s D — 
ev e/ Pp; 4 e/o a 











8 R cos 6 cos 6 dé “ (p—R cos 8) dé "Rp sin? 6 dé 
18s — { ie Di Be get. 
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Integrating by parts, the last integral in the right member 
becomes equal to the first and we get 


8 R cos 6 pot cos 6 

IR q i d= nf —— dé. 

é " R cos 8 do= - (sz—{) cos 6 dé 
dz 4 D i Ds : 


This . us 


= a te | ie { (p— xe) EP? 64 dk , 
eo es | ie 7 —_ ny £05 © aOdpdt , 


and we can therefore write 





Further 














cos w 8 

az= R & 
sin wo 80 

Hy=——- & (3) 
cos w 8® 

me oe 


From this it is easy to find the lines of magnetic force in space 
outside the whirl. These are the integral curves of the system 


dx _ dy _ dz 


H, H, d,’ 
where H,, Hy, and H, are the components of the magnetic force, 


parallel to the X-, Y-, and Z-axes respectively. 
But we have 


oer y 
Hr=> Het% Hy 
x 


Hs=5 H,—%, Z. 
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and if we take R and ¢ as variables instead of x and y, we have 


xdx+ydy=RdR , 
axdy— ydx= R'd¢ . 
Hence 
dx _dy_ xdx+ydy sa ady—ydx _ dz 
H, Hy xH,+yH, xH,—yH, 4H; 





and from the last three equations we get 


which are the differential equations of the lines of force in the vari- 
ables R, @, and z. 
Now if we substitute in the equation 


dR _ Rae 
Hr Hz 


the above-developed expressions for Hz and Hy, we get 


oF = —cot w dd, 


which, integrated, gives us 


R= Be-o«, (4) 


B being a constant of integration. 
On the other hand, if we substitute in the equations 


aR _ds 
He A, 


the values of Hr and H,, we get 
bm 5 
IR dR+> dz=o 
and ® being a function of R and z only, we can integrate and 
obtain : 
®=C (5) 


where C is a constant. 





Sa aE pe 


Sa ————— 


—~ 








— = 


6 aS SR ae a eS aS te 


a ee ee SS eee 


a 
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From equations (4) and (5) we have the following fundamental 


theorem: 
The lines of magnetic force due to the spiral whirl are curves in 


space lying on the surfaces of revolution ®=C. Their projections 
on the plane of the whirl are logarithmic spirals cutting the current- 
lines of the whirl everywhere in right angles. 


IV. FIRST APPLICATION TO THE THEORY OF SOLAR 
VORTICES 


Let us first consider the various cases that can occur for differ- 
ent signs of w,v, and e. In order to find the signs of the integrals, 


let us consider 
" cos 6.40 
. aes 
° 


D?= R?—2Rp cos 6+p?+(z—£)?. 


where 


This integral may be written 


" cos 6 “ cos 6 d@ " cos 6 dO 
{ Di w= | aaa | —_— 





By taking 0,=a—8@ as variable in the last integral, we get 


“008 iy_ { 200s 0d0_ ( 3 cos 6, dé, 
° Ds ° ° Ds ° D; 


D:= R?+2Rp cos 9,+p?+(s—{)? . 





where 


But for @=86,, and different from 2/2, we have D,>D, and there- 
fore the first integral will be greater than the second, that is, 


" cos 6 dé 
{ — PO 

















RESEARCHES ON SOLAR VORTICES 361 


On the other hand, ¢ and v have the same sign throughout the 
whirl, and if we are on the positive side of the X Y-plane and above 
the whirl, z—¢ will be positive. Farther along the Z-axis we 


have 
Pa eee Sd an ee 
3X10” 


Hence, from this equation and from equation (1) we conclude 
that: 

Her and H,, will have the opposite sign to ev, and Hy the same 
sign as ew?. 

We will apply this to solar vortices around sun-spots, assuming 
that the hydrogen flocculi represent lines of magnetic force due 
to a whirl in a lower level of the sun’s atmosphere. We will choose 
the Z-axis along a solar radius, with its positive direction pointing 
away from the sun’s center and toward the observer. For differ- 
ent signs of ©, w, and v, we have the 8 combinations shown in 
Table I, with corresponding signs for the components Hp, H,, 
and Hp. 























TABLE I 
TIE ee Ver eee + + ~ + ow 
Bathe oS NE eee aa sie Sale - + 
PAGES APSE EEN EA | + - de lt ~ - 
~ LAIR BS SSE | + - - | + + + 
Recess coke bau ankles oe i =—. | + + | + - + 
Yara bop. Gee oa ekee sao. — + _ - + - 
ERIC SLT Lom | 4 + 4. + 
Apppnapiets some eee eee ee ee ca | 











This gives the diagrams shown in Fig. 8, where the whirl is 
seen from a distant point outside the sun. The current lines are 
plain curves with arrows in the direction of motion, and the pro- 
jection of the lines of force situated between the whirl and the 
observer are dotted lines with arrows in the direction of the mag- 
netic force. When these last arrows point to the center, H,, is 
negative and the magnetic force along the Z-axis is directed away 


oa 
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from the observer; when they point out, that force is directed 
toward the observer. 


eav ft 








+++ * 
/ 
A | 
a | 
2. 
toa * 





- 
— 
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We now have to find which of these 8 cases corresponds with 
the observed facts. This can be done by combining Hale’s classical 
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researches on the Zeeman effect with the recent investigations on 
radial motion in sun-spots by St. John." 

In his paper ‘‘The Probable Existence of a Magnetic Field in 
Sun-Spots,”’ Hale states the following facts: “The magnetic force 
in the center of the whirl is directed toward the observer when the 
direction of motion indicated by the hydrogen flocculi whirl is 
clockwise, and directed away from the observer when the direction 
of motion is counter-clockwise.’” 

The meaning of the definition of clockwise and counter-clockwise 
motion will be seen’ in Fig. 9. There is here no question about 
direction of motion along the current- 
lines; the definition is purely geometri- 
cal. According to this, cases Nos. 1, 4, 
6, and 7 must be rejected: On the 
other hand, through the recent re- 
searches of St. John,‘ the existence of a 
whirl at a lower level giving the Zeeman fy 
effect and causing the observed arrange- ee 
ment of the hydrogen flocculi at the Clockwise 
higher levels is rendered extremely 
probable. In this low-level whirl, 
probably situated at or below the levels 
oo and ooo in the figure in St. John’s 
paper, St. John finds a motion outward 
from the axis of the whirl. This ex- 
cludes cases Nos. 2 and 3, and thus 
only cases Nos. 5 and 8 remain as cor- = 
responding to observed facts. In both Counter-clockwise 
cases € is negative, that is to say, we Fic. 9 
come to the same conclusion as Profes- 
sor Hale: The Zeeman effect is due to a whirl of negatively charged 
electric particles. 

This whirl, however, is distinct from the observed hydrogen 
whirl, and is situated in a much lower level of the sun’s atmosphere. 








t Mt. Wilson Contr., No. 69; Astrophysical Journal, 37, 322, 1913. 
2 Mt. Wilson Contr., No. 30; Astrophysical Journal, 28, 315, 1908. 
3 Op. cit. 4 Mt. Wilson Contr., No. 69; Astrophysical Journal, 37, 322, 1913. 
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The sign of electricity in the whirl being determined, we have 
to choose the probable value of the angle w, which determines the 
form of the current-lines. This can be done by comparing Fig. 10 
with a series of photographs of hydrogen flocculi whirls taken at 





-$0° 





FIG. 10 


Mount Wilson. It will thus be seen that values of w lying between 
15° and 25°, or between —15° and — 25°, agree well with the reality; 
the dotted lines must then be identified with the hydrogen flocculi. 
The current-lines and projection of lines of magnetic force (hydrogen 
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fiocculi) for the values w= 20° and w= —20° are shown on a larger 
scale in Fig. 11. This agrees well with the strong magnetic effect 
shown by solar vortices where the structure of the hydrogen flocculi 
is nearly radial; for the more nearly radial this structure is, the 
more nearly circular are the current-lines of the whirl. 

In order to get some idea regarding the amount of electricity - 
in motion in the whirl, more quantitative calculations are necessary. 


wa — 20° w= 20° 





——— Gwent lines of the whirl om the lower levele. 
--s---- Prozection of the lines of magnelie force [hydwogen. floceuli/) 
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For the sake of simplicity we will put forward another hypothesis con- 
cerning e and v as functions of p, as will be seen in the next paragraph. 


V. HYPOTHESIS REGARDING THE PRODUCT €p? AND THE CORRE- 
SPONDING FORMULAE FOR NUMERICAL CALCULATIONS 


We will now, for the sake of simplicity, make the following 
assumptions: 

1. The density of the gas in the whirl shall be constant through- 
out the whirl, and the charge ¢ proportional to the density; that 
is, € is supposed to be a constant. 

2. The velocity v shall be inversely proportional to the distance 
p from the axis. This hypothesis, as already known, was made by 
Guldberg and Mohn for the outer part of terrestrial cyclones. 
Thus epv will be constant throughout the whirl. Let us put 


Epv me 
3X10" 
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Hence 
* (" ('* (2-0 cos 6 
Hr=—2¢ cos -{ if J aah Nn d0dpdé 
D3 
° PI °o 
(*h (*p, 
Hy= 2csinw " = ot ED : d6dpdé 
Jo ed, eo 
(*h (*p 
. 6— 
H,= 2¢ cos» * R cos eth d0dpdé 
Ds 
Jo e/P, eo 
and 
cos 6 
omar | : i jp ddpdé 
where ; 


D?= R?—2Rp é0s 6+p?+ (z—)?. 


We now have 


: dt=7,+ constant ' 


fe 
R cos 6—p I 
—— dp= pt constant ° 





ed) 


Hence we may write 


Hr=-—2c¢ cos w [f(h)—f (0)] 
Hy= 2c sin o[f(h)—f(0)] 
H,= 2c cos [g(p.)—g(p:)] 


where 


p 7 h 7 
. 6 d6d 
jo~ | { Abdo d= { 7 


(6) 
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But we have further, /m being the natural logarithm, with base 
e=2.71828 ....: 


fe 








oP In (p—R cos 6+D)+ constant , 
, 
re dt 
5” —In (s—f+D)+ constant , 
y, 
which gives 
f@) =f(p2, g) —f(p:; é) ’ 
g(p)=g(p, h)—glp, 0), 
where 
(nr 
f(e,0)= In (p—R cos 6+D) cos 6 dé 
e/Jo (8) 
(*n 
g(p,f)=— In (2—{+D) dé 
e/Jo 


Substituting into the expressions for the components, we 
finally get 


Hr=—2¢ cos » [f(p,, h)—f (p1, h)—f (p2, 0) +f (p:, 0)] 
Hy= 2c sin » [f(p., h)—f(p:, h)—f(p2, 0) +f (m1, ©)] (9) 
H,= 2¢ cos » [g(p., h)—g(p:, h)—g(p., 0) + (¢:, ©)] ) 


For each system of values, R,2,9,f, the integrals f(p,¢) and 
g(p,¢) can now be easily calculated numerically. Integraphs and 
planimeters may be employed, or the approximate formula of 
Simpson, which is well known. 

We thus have the means of calculating in every case the mag- 
netic field in space outside the whirl. 

If the point (x,y,z) be on the axis of the whirl, R is zero and 


Di=r+(s-0, 
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which gives 
t (p,g)=ln to) | cos 6 d6=o0, 


g(p,6)=—m In (2—£+D). 
Hence Hr=o, Hy=0, and H,=H., where 


(2—h+V (2—h)*+p2)(s+1 z+pi) 


—— (10) 
(s—h+V (s—h)?+p3)(s+V 27+ 3) 


Ay»=—27¢ coswln 


It will be interesting to obtain a verification. Let h and k 
be infinitesimal, and 


=ptk, pi=p, w=o. 





If we put oe 
V pP+e=u 

we have ie 

s—h+V (—hy+ (pt hP= 

pies Sal id... ps 
(e+w)| 3 5+ gta Maoh + (s+u)us ae 

Teta a 4 

Hence 





s— —h+1 (s—h) +( pth)? _ p p u u?-- us 2 
thei agent } ———— aan 
er h+V (c—hy+p: (a+u)u a a 


This gives 





In 2. Verp  -» p—e ue t2 
stVe+(pth) — (etu)u  e+u)w 





i ee 


and finally 
H,=— 27 p hWk+.... 
u3 
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where 
epv 


6 3X10" ° 


If we suppose e and 2 to be positive, e#kv electrostatic units 
will pass through the cross-section of the current in one second, 
which corresponds to a current of 





. 
~ 3X10? 
amperes. If we introduce this value, 
. 2m? 
H,= oan 4 4 4. oe 
10 “8 


and this first term represents exactly the magnetic action in the 
point (0,0,z) of the circular current whose section is the infinitesimal 


Z 








bj 
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product Ak, and where the direction of the current is that of 
increasing @ (see Fig. 12). 
VI. SPECIAL CASE IN WHICH THE THICKNESS OF THE WHIRL IS VERY 
SMALL COMPARED WITH ITS DIAMETER 
We will now consider the special case in which the thickness 
of the whirl is very small compared with its diameter. Then for 
points (x,y,z) not very near the whirl, we can put 


f(p,h) =f (p,0) =hf"(p,0) ’ 
g (p,h)—g (p,0) =hg’ (p,0) , 
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where f’(p,o) and g’(p,o) are the partial derivatives of and 8 for 
6f 6¢ 
the value (=o. 
Now 


’ Pe et. z—¢ cos 6 dé 7a) °* 46 
f (p,f) = { p—R cos 6+D - eae Th » & i= | D 


and 


: r cos 6 d6 ” (p—R cos 8) cos 6 dé 
f'(po)=—2 peepee —— mg lt Bh ti! By Be 
, 4 (ep—R cos 6+-D,)D, a (R? sin? 6+27)D, 


g ood f . ; 


where 
D2= R?—2Rp cos 6+-p?+2?. 


This must be substituted in the formulae 


Hr=—2ch cos » [f'(p,, 0)—f"(p:, ©)] 
Hy= 2ch sin » [f'(p,, o)—f"(p:, 0)] (11) 
H,= 2ch cos w [g’(p,, 0)—g"(p:, ©)] 


The same result can also be obtained by starting from formula 
(6), and considering / as infinitesimal. Then 


p wT 
’ 6 
Hr=—2ch cone [ j — d6dp 
Py ° 


and an integration with respect to p gives the foregoing result. 
The function @ can easily be found in this case. We have 


P2 7 
cos 6 
i= 6, 
2ch ef j Dz d6dp, 


fg =Iln (p—R cos 6+-D,)+ constant, 





and as 
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we get 


= 2chR [ f(p2, 0) —f(p:, ©)] (x2) 


where the function /(p,f) is given by equation (8), when D=D,. 

The integrals occurring in the formulae for the components 
can be expressed by elliptic integrals of the first and second kind; 
but only the reduction of the integral for g’(p,o) is sufficiently 
simple to be of any importance in the present problem. The 
integral in this case simplifies by the well-known substitution’ 


as follows: 


. d6 2 {; dy (a3) 
ee en I 
* R?—2Rp cos 6+p?+2 A ce 1—k? sin? w 3 








where 





Page A=V (R+p)*+2". 

As this integral is well tabulated in Legendre’s Traité des 
Fonctions Elliptiques (Paris, 1826), 2, 223, the calculation is rather 
easy. 

The other integrals are most conveniently computed by Simp- 
son’s formula, or by machines. 


VII. NUMERICAL APPLICATION TO A SOLAR VORTEX 


We will now apply the preceding formula to an ideal solar 
vortex, in order to see how many of the observed facts can be ac- 
counted for. 

According to the researches mentioned above by St. John, the 
vortex giving the Zeeman effect is probably situated at a level 
of the sun’s atmosphere equal to or lower than the level Fe, 00, 
where the pressure is probably about 10 atmospheres. As the 
pressure increases very rapidly in descending into the sun’s atmos- 
phere, it seems probable that the thickness of the outflowing whirl 


* See, for instance, Appell, Traité de Mécanique Rationnelle, 3, chap. xxxv. 











372 CARL STORMER 


in question can be only a very small fraction of the diameter indi- 
cated by the hydrogen flocculi whirl in the upper levels. This is 
also analogous to terrestrial cyclones, whose thickness is very small 
compared with their diameter. In other words, # must be very 
small compared with p, As to the size of p, compared with h 
and p2, we can make a probable assumption by starting from the 
fact observed by Hale, namely, that the magnetic effect of the 
vortex decreases very rapidly upward along: the axis of the vortex, 
so that there is only a very small fraction of its original value 
when we come to the level of 5000 miles above the photosphere. 
In order to do this, we will start from the formula 


(s—h+V (—hy +p) (2+ VF +p) 











Aw~=—27coswln SEY : F 
(s—h+V (s—h)?+p)(2+V 27+ ?) 
This gives 
dH x ( I I I I ) 
—~=2"¢Cosw 4 ————— — ——— =—- ———— > 
ds (Va+pi V(s—h)*+p2 Ve+p? V (s—h)?+97) 


The geometrical signification of the radicals will be seen in Fig. 13. 
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On studying this figure, it will be seen that if we choose h 
small compared with the elevation of the high levels of the sun’s 
atmosphere, and, on the other hand, suppose p,; to be of the same 
size as, or smaller than, 4, the magnetic force at these high levels 
will be only a very small fraction of what it is at the level z=/. 
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We will choose as a numerical example 
h=p,=500 km=5X10’ cm, p,=50,000 km=5 X10? cm. 


The last value is chosen to correspond with the size of the observed 
hydrogen whirls around sun-spots. 

Further, we take w=20°. .With regard to the constant c, we 
will start from the value of about 3000 gausses observed in sun- 
spots by the Zeeman effect, and suppose that this is the value of 
H,. for z=h. Wethen get c=—583.1. Starting from these values 
we find the following values of H,, along the Z-axis: 


sin km H,, in Gausses sin km H,.. in Gausses 
MG 6 dab aeatan ti enaad 3000 | RRR MRT, FEL! 99> 142.6 
GS: § 3 oaks Sewiwhedkie IQOI CA. D ANS S000 cb eee 55-33 
Pa é-s ukesevencaevenee 916.7 Sis 6 in.eeee teks we 28.46 
RE Se ee 583.2 Ge oss cane nied 16.48 
ere e 318.4 | ee er eee 10.62 


A probable system of values for h, p;, p2, w, and c being found, it 
is of the greatest interest to compute the magnetic field above the 
spot. We have done this for each point encircled in Fig. 14. 
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We have used the formula given in Section VI for the case in which 
h is small compared with p,. In the following paragraphs we will 
give these computations in detail. 


VIII. COMPUTATION OF THE COMPONENTS H; AND He 
We have 


H,=4he cos w [2-<] 


SRE SEES es 
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where the subscripts 1 and 2 denote the values for p=p, and p=p, 
of the functions 


T 








ere a 
A=V(Rteyte, T=) 7 sing’ 
where 
wert eee 
N(R+e) +3" 
and further 
h =5X107 w= 20° 
Pr=5 X10! c= —583.1. 
P2= 5X 10° 


By means of Legendre’s well-known tables, we find the results 


given in Table II. 
TABLE II 


VALUES oF H: 


Rin Km 








| 
zin Km —= ———— ‘me ee soe — a 

| 10,000 | 20,000 30,000 | 40,000 | 50,000 | 60,000 | 70,000 | 80,000 | g0,000 | 100,000 
50,000... .| 9.47| 7.89) 5.84) 3.84 |4+2.22/+1.09/+0.39/+0.01,—0.17/—0. 25 
40,000..../ 14.98) 11.74) 8.01) 4.64. |-+2.18|+0.68|—0.08)—0. 39 —0.48,—0. 38 
30,000. ...| 24.86) 18.09) 10.94) 5.22 |+1.51|/—0.31, —0.92|—0.99|—0.88 —o. 73 
20,000.... 44.84, 28.27) 14.49) 4.85 |—0.65|—2.34 bape Wha. —1.31|—0.99 
10,000... .| 87.65! 42.00) 17.72) 2.33 |—6.58)/—5.95 —3-76|—2.43|—1.66|—1.18 





The calculation of the other component, Hr, has been much 
more arduous. We have used Simpson’s formula for the approxi- 
mate computation of the definite integrals: 


b 
: h 
( f(x)dx= : [Yo 4 yt 2V2t4yst2Vgt «~~ H2Vn—2b4¥n—r In] » 


a 


where the interval of integration b—a is divided into m parts, 
each of them equal to /, and where yo, ¥:, V2 - - - - Va—1) Vn are the 
values of f(x) for the points of division. Now 


Hr=—2chz cos » [Y,—Y;] 
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where Y, and Y, are the values for p=p, and p=p, of the integral 


r-{ (p—R cos 8) cos o 0 





(R? sin? 0+-2*) D ‘ 


and 
D?= R?—2Rp cos 0+ p?+2. 


Here the interval from o to x has been divided into 12 parts, 
which has involved the calculation of the function under the sign 
of integration for each combination of the values mentioned above 
of R, z, and p, with the values 


0=0°, 15°, 30°, 45°) . - . » 150°, 165°, 180°. 


As the value of Hz is not affected if we divide each of the quanti- 
ties R, Z, p:, p2, and h by the same number, that is, if we choose 
another unit of length, we have chosen / as this unit. The values 
of the quantities will then be 


h=1, pr=1I, p2=I100 
R= 20, 40, 60, 80, 100, 120, 140, 160, 180, 200 
Z= 20, 40, 60, 80, 100. 


The calculation has been made with the “‘ Brunswiga”’ machine. 
The results for the function 


(e—R cos 4) cos 6 
(R? sin?0+27) D 
were integrated by Simpson’s formula. Multiplication by the 
constant — 2chz cos w gave finally the results shown in Table III. 


TABLE III 
VALUES OF Hr 











sin Km wee neat ier 
| 10,000 | 20,000 | 30,000 | 40,000 | 50,000 60,000 | 70,000 | 80,000 | 90,000 | 100,000 
50,000. 2.12 | 3.69 | 4.44 | 4.44 | 3.95 | 3.20 | 2.47 1.37 | 0.73 





40,000....| 3.82 | 6.37 | 7.24 6.79 | 5.61 | 4-20 | 2.99 

-+| 7.50 |11.65 |12.25 10.57 | 7.93 | 5.30 | 3.36 
20,000... .|17.04 |22.71 |21.17 |16.58 |10.92 | 6.12 3.33 
10,000... .|49.68 \45 .86 35.56 25.48 |13.83 | 5.30 | 2.28 


86 

08 | 1.45 | 1.02 
.15 | 1.40 | 0.91 
gO | 1.15 | 0.73 
16 | 0.66 | 0.42 














“x -# YH 











ESSE 


a 


acne ee eget 


es a oe eee 
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By the formula 
H,=—Hetano, 


where 


tan w=tan 20°=0. 36397, 


we find the tabulated values of H, (Table IV). 


TABLE IV 


VALUES OF Ho 

















Rin Km 

s in Km | ~ ‘RES SERA ERE Benen 

10,000 20,000 30,000 | 40,000 50,000 | 60,000 70,000 | 80,000 | 90,000 | 100,000 
50,000 .|— 0.77/— 1.34/— 1.62|—1.62|—1.44|—1.17/—0.90 —0.68|—0. 50 —0.27 
40,000 .|— 1.39/— 2.32/— 2.64|/—2.47|/—2.04/—1.53|—1.09|—1.76|—0. 53) —0.37 
30,000 .|—  2.73|— 4.24/— 4.46|/—3.85|—2.89|—1.93|/—1.22|—0. 78/0. 51, —0. 33 
20,000 .|— 6.20/— 8.27/— 7.70|—6.03|/—3.97|/—2.23/—1.21|—0.69] —0.42| —0. 27 
10,000 .| — 18.08} — 16.69 — 12.94] —9.27|—5.03| —1.93|—0.83 —0.42/—0.24)—0.15 


By means of these values we can calculate the magnetic force 
H=V Hyt+H:+4j. 


TABLE V 


VALUES OF H 


























Rin Km 
zin Km ————— = . ] - 
| 10,000 20,000 30,000 | 40,000 50,000 60,000 | 70,000 $0,000 90,000 | 100,000 
} | | 
50,000....| 9.74) 8.80) 7.51) 6.09) 4 75| 3.58 | 2.66 1.98 | 1.47 | 0.81 
40,000....| 15.52) 13.56) 11.11) 8.59) 6.36) 4.52 | 3.18 | 2.25 | 1.61 | 1.15 
30,000....| 26.10) 21.93} 17.02] 12.40} 8.57] 5.65 | 3.69 | 2.49 | 1.73 | 1.21 
20,000....| 48.37) 37.20| 26.79] 18.30] 11.64] 6.92 | 4.18 | 2.66 | 1.79 | 1.26 
10,000... .|102.4 04-39) 41.78] 27.21| 16.13] 8.20 | 4.48 | 2.73 | 1.81 | 1.26 





The component in the meridian plane through the axis of the 
whirl has the components Hp and H, and is everywhere tangent 
to the curves ®=constant, where ® is the function given by 
equation (12). 
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The component in the meridian plane for some of the points 
(R,z) in the tables mentioned above is shown in Fig. 15. It will 
be seen that the magnetic force decreases rapidly upward from ° 
the whirl. 

Seale of magnetic force 


we ces. 


. ‘ \ \ I / / / 

o . \ \ \ 1g / / / - m 
o ~ + \ \ ~ a / vA a - = 
la, Mags / fee a 
fooge han 





50008 ham 


Magnetic fil of solat vortex. 
Fic. 15 


IX. CALCULATION AND CONSTRUCTION OF THE LINES OF FORCE 


We have seen that the lines of force in space are the curves 
of intersection between cylinders through the logarithmic spirals 


R= Be-¢ te 
and the surfaces of revolution, whose meridian curves are 
=C. 


In the case in question we have 


= 2 chR [f (p., o)—f (p:, 0)] 


where 


riosd= { In (p—R cos 6+D) cos 6 dé, 


D?= R?—2Rp cos 6+)p?+2'. 
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In order to draw the lines =constant, the following method 
has been used: For the computation of the integral by Simpson’s 
’ formula, we have first computed the function 


log (ep—R cos 8+D) cos 8, 


where the logarithm is to the base 10 and where p, R, z have the 
values given before, with / as unit of length. 
If we denote by 


Veo ar Pas 2 0 0 o Yes 
the values of this function for 
O6=0, 15°, 30°, 45°, - . . - 165° and 180° 
and S, and S, the values of the expression 
Vor 4Vit2VzF4yst +--+ F4¥utyu 


for p=p, and p=p.,, we then have approximately 





(3 
= S,= log (p:x—R cos 0+D,) cos 6 dé 
3 ope 
(? x 
- 5,= log (p:—R cos 6+-D,) cos 6 dé 
eJo 


where D, and D, are the values of D for p=p, and p=p,. As 
log A=log e - In A=0.43429/n A 


where ¢ is the base of the hyperbolic logarithms, we get approxi- 


mately 
ach 
Oii<en— © Zt) 
18 log e R (S.—S1) 

Thus the function R(S,—S,) is a function of R and of z, which 
‘ differs from the function ® only by a constant factor and this even 
if we express R, €, p:, p2, and /# with / as unit of length. 

Thus the curves @=constant are identical with the curves 
@,= constant, where &,= R(S,—S,), if Simpson’s formula is assumed 
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to be accurate enough for our purpose. From this we have then 
drawn the curves with R as abscissa and ®, ordinate, and where z 
has a constant value for each curve. 

The points (Ro, 2) and (Ro, 2) where ® has a given value, 
®,, can then be found graphically by the intersection of the curve 
and the line =, (Fig. 16). Similar curves have been drawn 
with z as abscissa and ® as ordinate, for controlling the results. 
This has finally given enough points for the construction of the 
curves ®,= constant. 


, 




















R, R, R 


r 2 


Fic. 16 


In order to control the results, the function ®, has also been 
calculated for z=o, in which case the integrals can be reduced to 
elliptic integrals of the first and second kind by a partial integration. 

The curves ®,= constant are seen in Fig. 17, and the directions 
of the components of magnetic force in the meridian plane are also 
given. It will be seen that the agreement is fairly good, in spite 
of the approximate method of computation and construction. 

With the aid of the curves ®,=constant and the logarithmic 
spirals 

R= Be? ete , 


where B is a constant and w=+20°, I have constructed wire 
models. of the lines of magnetic force in space over the whirl. The 
dimensions are the same as for the case already mentioned. 
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PLATE VII 


STEREOSCOPIC VIEWS OF MODEL SHOWING DISTRIBUTION OF LINES OF MAGNETIC FORCE FOR w= — 20 
CuRRENT LINES IN LOWER LEVELS OF WHIRL ARE TRACED ON THE BASE OF THE MODEL 
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In Plate VII are seen some stereoscopic representations of these 
models for the case w= — 20°, viewed from various points. On the 
base are drawn the current-lines, and the lines of force in space are 
supported by dark pins, standing perpendicular to the plane of the 
whirl. 

If we compare the curves ®,=constant with the well-known 
spectroheliograph pictures of calcium prominences around a sun- 
spot at the sun’s limb, shown by Professor Slocum in his paper* 
“Attractions of Sun-Spots for Prominences,” a very striking 
resemblance appears. It seems as if these prominences are arranged 
along lines of magnetic force from a whirl round the sun-spot. 


X. THE MAGNETIC FIELD ROUND A BIPOLAR SPOT-GROUP 


Let us now consider two whirls and their combined magnetic 
field of force. For simplicity we will suppose that each of the 
whirls is of the kind just mentioned, and that they are lying in the 
same plane with parallel axes perpendicular to the plane. The 


a 








' Fic. 18 
Se 
Fic. 19 


magnetic resultant force in each point of space over the whirls can 
then be found by the parallelogram law, because we know the 
force due to each of the whirls considered separately. 

The graphical method used for finding the resultant components 
at a given level above the plane of the whirls is as follows: on a 
sheet of tracing paper are marked the components of the magnetic 


* Astrophysical Journal, 36, 265, 1912. 
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force due to the whirl A (Fig. 18), and on another sheet the com- 
ponents due to the whirl B (Fig. 19), both components for points 
along a radius vector in the selected level z=constant. The two 
sheets are then placed with the points A and B coinciding with the 
centers of the whirls (Fig. 20), and turned to produce coincidence 
of the points M and N on the sheets. On a third underlying sheet 
of tracing paper the points M and R are marked by pricking with a 
pin, and the resultant force MR is then constructed. 











Fic. 20 
X 
M 
P 4 
a i 
A r “B Y 
FIG. 21 


With regard to the calculation of the resultant force the neces- 
sary formula may be developed as follows: Let us consider a given 
plane parallel to the plane of the whirls, and let A and B be the 
points where the axes of the whirls intersect this plane. Let 7 be 
their mutual distance. Let M be a point in the plane and Hg, 
H,, and H, the components at M due to the whirl with axis A, and 
Kr, Ky, and K, the components due to the whirl B. Let us finally 
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denote the distances AM and BM by p and g and the angles MAB 
and MBA byaand f. Then 


cos on tO , 
2pr 
cos pa tr—P : 
2gr 
If we choose a system of rectangular co-ordinates with the 
origin in the center of the whirl A, and the X-, Y-, and Z-axes 


parallel to AX (perpendicular to AB, see Fig. 21), to AY, and to 








Fic. 22 


the axis of the whirl, and if we denote by P,, P,, and P, the com- 
ponents of the resulting magnetic force at the point M, parallel 
with the X-, Y-, and Z-axes, we find 


P,= Hp sin a— Hy cos a+ Kp sin B+Ky cos B, 
P,y=Hrcosa+Hy sin a—Kp cos B+K;y sin B, 
P,=H,+K;. 
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For the practical computation of P,;, P,, and P,, sin a, cos a, 
sin 8, and cos 6 can first be calculated. If, for instance, r= 10,000 
km, the distances p and g have the values: 


10,000, 20, 


000, . . 


. « 90,000 and 100,000 km, 


and we obtain the values in the following tables corresponding to 


the numbered points in Fig. 22. 


The values of these functions will 


be useful for the calculation of P,, Py, and P, in special cases. 










































































TABLE VI 
g = 100,000 
p 
l l ; | ak a ‘eve: 
100,000} 90,000 | 80,000 | 70,000 | 60,000 | $0,000 | 40,000 | 30,000 | 20,000 | 10,000 
ee ©.8660]0.8930/0.9165/0.9367 ©.9539|0.9682/0.9798/0.9887/0..9950/0.9987 
ae ©. 5000/0. 4500/0. 4000/0 . 3500/0. 3000/0. 2500/0. 2000/0. 1500/0. 1000/0 .0500 
| Se ©.8660)0 . 8037/0. 7332/0.6557/0. 5723/0. 4841/0. 3919/0. 2965 0. 1990/0.0998 
| ©. $0000. 5950|0. 6800/0. 7550/0. 8200/0. 8750}0.9200]0.9550/0.9600)0.9950 
Point no...|  t 2 | 3 | 4 | 5 6 | 7 8 9 10 
TABLE VII 
§ =90,000 
- p 
90,000 | 80,000 70,000 60,000%| 50,000 | 40,000 | 30,000 20,000 
ee 0.8315) 0.8549) °.8741| 0. 8888| 0.8980) 0.8992) 0.8843) 0.8181 
ae 0.55506) 0. 5187) 0.4857) 0.4583) 0.4400] 0.4375) 0.4668) 0.5750 
Ns aid omen 0.8315) 0.7599} 0.6798) 0.5924) 0.4988) 0.30995) 0.2947) 0.1814 
See 0.5556) 0.6501) 0.7334) 0.8056) 0.8667) 0.9167) 0.9556) 0.9834 
Point no... .. —) | 13 | 14 | 15 | 16 17 18 
TABLE VIII 
g = 80,000 
; ———————————— 
| ab a = 
80,000 70,000 60,000 50,000 40,000 30,000 
igo Gta che ecasshanid hk od 0.7806 | 0.7945 | 0.8000 | 0.7924 | 0.7599 | 0.6613 
DMbie.4 Perdeveeh er ees 0.6250 | 0.6072 0.6000 | 0.6100 | 0.6500 | 0.7501 
ah: and dia Sinden ees | 0.7806 | 0.6953 | 0.6000 | 0.4953 | 0.3799 | 0.2480 
I ig Sinsin wah gtotal sen ot | 0.6250 | 0.7187 | 0.8000 | 0.8687 | 0.9250 | 0.9687 
SP Ms. i ceo y at’eia » | 19 20 21 | 22 23 24 
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TABLE IX 
§ = 70,000 
? 
70,000 60,000 50,000 40,000 
Gs ep aiean ans Gaede shee bod 0.6998 | 0.6887 0.6500 0.5464 
Ms obs as ae CREE eens | 0.7143 | 0.7250 | 0.7600 0.8375 
GE 65 $4.0 6k vsasabuaton 0.6998 | 0.5903 | 0.4642 ©.3122 
GUD kesh od 6ensteneeaexs | 0.7163 | 0.8072 | 0.8857 ©.9500 
Pets ca ewdesassena 25 26 27 | 28 
TABLE X 
g = 60,000 
ES — a —_ — — 
p 
WepeoeneNeeMaGS = 
60,000 50,000 
iii wis doetniwes | 0.5528 0.4560 
Got ad cvauceus 0.8333 | ©. 8900 
GE ok bau oa tee 0.5528 ©.3799 
OEE es 0.8333 0.9250 
PUES sib tcxedas 29 30 





XI. SOME SIMPLE CASES OF BIPOLAR SPOT-GROUPS 


There are two interesting cases which we will study in greater 
detail, namely, the case of two identical whirls, and the case of two 
equal and opposite whirls. In both cases there is a certain sym- 
metry in the arrangement of the lines of force, and their projection 
on the plane of the whirls. 

First case: two identical whirls.—Let us consider the compo- 
nents P, and P, and their resultant P, in a plane parallel to the 
plane of the whirls. Let C (Fig. 23) be the middle point of the 
line AB joining points A and B where the axes of the whirls cut the 
plane. Let M and M’ be two points in the plane situated sym- 
metrically with regard to C, and let us denote by accents the 
letters corresponding to the point M’. We then have 


a=—B , ‘=—a 
Kin . Sint 
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which gives 
P!l=—Hpsina+H, cos a—Kpr sin B— Ky cos B=—P, 
P)\=—Hpr cosa—Hy sina+Kr cos B—Ky sin B=—Py, 
i.e., the component Pj, in the point M’ is obtained by turning the 
component P, in the point M through an angle of 180° around C as 
center. 

















A jc B 
Fic. 24 


In Fig. 25 are seen the directions of the components P, in 
a plane, corresponding to the level 20,000 km over the two 
whirls in the sun’s atmosphere. For the two identical whirls 
we have chosen case No. 5 of Fig. 8, with 
h=500 km, w= 20° 
pPr=500 km , c=—581.3 
P2=50,000km, r=100,000 km. 
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Second case: two whirls corresponding to cases Nos. 5 and 8 in 
Fig. 8.—For the whirl A we have ¢ negative and w and 2 positive, 
and for the whirl B, 


, , 


=e, w=—w, v=—9, 


consequently oo 
c’=pe'y’=—c. 


Further, for both whirls 4, p,, and p, are supposed to be the same. 
Let us consider (Fig. 24) as before a given level above the plane 
of the two whirls, and let CD be the perpendiculfr through the 
central point C of the line AB. Let M and M’ be two points 
situated symmetrically with respect to this line, and let us denote 
by accents the letters corresponding to the point M’. 
We have 


a’=B, f'=a, 
Further by formula (11), Hp positive, Hy, negative and 
Hp=—Kr, Kr=—Hr 
Hj= K,, Kj= Hy 
we find 
P'=+Hz} sin a’—H} cos a’+ Kk sin B’+Kj cos p’ 
=—Hrsin a +H, cosa —Kpz sin 8B —Ky cos 8 =—P, 
and 


P} =H cos a’+Hj, sin a’— Kz cos B’+ Kj sin p’ 
= Hr cosa +H, sina —Kr cos 8 +Ky sin B =P,. 


Hence the straight lines passing through M and M’ in the direc- 
tion of the force lie symmetrically with respect to the line CD. 
If M is on the line CD, we have 


Hr=—Kr, Hy=Kz, a=B 


and 
P,=0, Py=2Rcosa+2Hy sina, 


i.e., the component P, is parallel to the line AB. 
In Fig. 26 are seen the directions of the component P, in 
a plane corresponding to a level of 20,000 km, in the sun’s 
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‘atmosphere, where the hydrogen flocculi are probably to be found. 
The dimensions and the distance of the whirls are chosen as in 
the first case. 

It is naturally of considerable interest to compare this with 
spectroheliograms of bipolar sun-spot groups. This will be done 
later, when more material has been obtained. 


XII. THE ELECTRIC CHARGE OF THE SOLAR VORTEX AND THE 
ELECTROSTATIC FIELD OVER IT. DIFFICULTIES IN EXPLAINING 
THE ABSENCE OF THE STARK EFFECT. HYPOTHESIS OF GAL- 
VANIC CURRENTS 


We have seen that the charge of the low-level whirl probably 
causing the Zeeman effect is negative. The value of the constant 
c will give us the means of determining this charge and of calcu- 
lating the electrostatic field in space over the whirl. We had, by 
definition (see Section V), 


c= PEP 
3X ro” ? 





where ¢ is the charge of electricity in 1 cubic cm of the whirl, 
measured in electrostatic units, and p and v two corresponding 
values of the distance from the axis and the velocity of a gas- 
particle of the whirl, measured in the c.g.s. system of units. 
€ was supposed constant throughout the whirl. 

If we measure directly the velocity at a given distance from the 
axis of the whirl, we can compute pv, and c being known from the 
Zeeman effect, we can find e. 

Let us make an estimate of the probable value of the charge 
e. As stated in his paper ‘Radial Motions in Sun-Spots,’” 
St. John has found, at the level Fe, oo and at a distance of 11,250 
km from the axis, an outward motion of about 1 km per second. 
If we assume the current-line to be a logarithmic spiral correspond- 
ing to w= 20°, the velocity will then be 


105 
v=-— 5 cm 
sin 20 





* Mt. Wilson Contr., No. 69; Astrophysical Journal, 37, 322, 1913- 
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i.e., about 3X105 cm. If we choose the value of c determined 
in Section VII, namely, c=—583.1 we find e=—0.052. From 
this we can compute the number WN of corpuscles to a cubic cm 
of the gas, as we know the elementary charge of a negative cor- 
puscle, which is*' —4.65X107*°.. We then find V=1.1X 108%. 


Z 


3 








Y 
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Let us now compute the electrostatic field above the whirl. 
Consider an element of the whirl bounded by the planes (Fig. 27) 


z=¢, z={+AZ, 
the planes . 

$=6, o=0+A0, 
and the cylinders 

R=p, R=p+Ap. 


The volume of this element is equal to pApA@Af¢ and its charge is 
epApAbds. 

Let 5£,, 6E,, and 5E, be the components parallel to the axes 
OX, OY, and OZ of the electrostatic action of this element on one 


*See Rutherford, Radioactive Substances and Their Radiations, § 20, Cambridge, 
1913. 
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electrostatic unit of positive electricity in the point (x,0,z) of the 


XZ-plane. Then 


x— pamet 


SE,= epApAdQZ , 


; 


~eeee 


dE,= epApAGQ , 


: 


bE, = * epApAddl . 
The action of an element lying symmetrically with respect to the 
XZ-plane is obtained by changing the sign of @. Here r, is the 
distance from the point (x,0,z) to the element. 

The components of the resultant action of the two elements 


will be 
25E, , Oo, 26E,. 


If we suppose the whirl to be of the same shape as in Section 
II, the components of its action on an arbitrary point (x,y,z) will 
be 


rP2 (°* (°" (R-» cos 6)pdpdtde 
Ep=2e a 
e/ Py eJo Jo 
(14) 
eeaoe(” (" (0 G_Sedeatas 
‘Sie Ds 











ep, e/0 eo 


Here Ez and E, are the components normal to and parallel to the 
Z-axis, lying in the meridian plane and with the same orientation 
as the components Hex and H, derived in Section II; further, 


D?= R?+-p?— 2Rp cos 6+ (s—{)?. 


If we introduce the potential 


of Sf 


sw sw 
8R’ 


we can write 


Er=-— E,=- 


&: * 
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The components of action of the whirl on an electric charge of 
e’ electrostatic units in the point (x,y,z) are consequently e’Epz 
and e’E£,. 

If the thickness / of the whirl is very small compared with its 
radius p2, we can use the formula corresponding to h infinitesimal, 


namely, 
"* (" pdpdé 
Waser | f - 
with 
B=, Bae, 
where 


D2 = R?+-p?— 2Rp cos 6+27. 


The expression for W can easily be transformed into another 
form that is better adapted to numerical computation. We have 


jz =D,+R cos 9ln(p—R cos 6+-D,)+constant, 


which gives 
W= W (p.) oat W (p:) ’ 


where 


° e (16) 
W(p)= wh | Dd0+ wie f cos 6 In(p—R cos 6+D,)d0 ) 


° ° 


Here the second integral is the same as in formula (12) for the 
function ®, and has been denoted by f(p,o). The first integral 
can be expressed by elliptic functions by the substitution quoted 
in Section VI, which gives 


| rae—a [ I—F sin? o dy : 


dX and & having the same definition as before. This integral 
can then be calculated by Legendre’s tables. 


—— = 


| 
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Let us now consider the components Erg and E,. We must 
first find the partial derivatives of W(p) with respect to R and 
to z. 

We have 


8W(p) _ "R—p cos 0 ry : 
— wh | ae d0+-2¢eh a cos 6 In(p Rcos 6+-D,)dé ° 


° ° 


But the partial derivative in the right member becomes 


. " R cos 0— (e+D.) cos? 8 
if cos 6 In(p—R cos oroater | | D.(p—R cos 6+D,) = 





° 


Integrated by parts, the first term of this becomes 


- j " (Ret+RD,) sin? 94g 





D.(p—R cos 8+-D,) 


° 


which, combined with the second term, gives for the partial deriv- 


ative 
"46 
° D 


SW (p) _ “cos 6 d0 
“_ inf Dp. ° 


° 


Consequently, 





If we use the transformation of Section VI, we further obtain 
this integral expressed in terms of elliptic functions in the form 


Er= 2eh{T (p.) ox T(p:) ] ’ 
where 


2p = dw Saas 4 b 
To=ie]e—¥ 1 = awe V 1—F? sin Za 


Finally we have to determine E,. By partial differentiation 


SW (p)_ "40 "cos 6 d0 
5 “=F 5T) Chast: 











° 
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Here we have the integrals which in Section VI we have denoted 
by /’(p,0) and g’(p,o). We find 

E,= — 2ehz[g'(p2, 0) —g (p:, 0)]+2ehR[f’(p., 0) —f"(p:. ©)] (18) 

By formulae (16), (17), and (18), the functions W, Ep, and E, 


can easily be computed. 
Along the axis of the whirl, 





W =ameh[V p2+2—V p?+2] 


Er=o (19) 











E,=2nch| 6 4 | 
V +p? V 2+)? 

Let us apply these to the solar vortex, studied in Section VII. 
We have 


h=5X10' cm, pPr=5X107 cm, P2=5X109 cm, 
€=—0.052 electrostatic units. 


Let us choose z= 10°, 2X 10°, 3X10, 4X 109, and 5X10 cm. We 
then find: 


Z= 10,000 km | 20,000 30,000 40,000 | 50,000 





| 
E;=—1.3X107 | — 107 | —o.8X 10! | —o.6X 107 —0o.5X 107 


Here, as we have seen, E, is the force acting on one electrostatic unit 
of positive electricity. 

The immense values of this force compared with the absence 
of any appreciable Stark effect over sun-spots' seem quite improb- 
able. Even with a much smaller convection whirl the observed 
magnitude of the Zeeman effect leads us, as Salet has already 
remarked,? to electric forces corresponding to several million 
volts per centimeter. We have tried to explain this paradox by 
assuming with Deslandres that the sun’s atmosphere is composed 
of different levels of opposite charges. If one of these levels is in 


* According to a letter from Professor Hale dated June 30, 1914. 
2 Bulletin de la Société Astronomique de France, 28, 220, 1914. 
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rapid motion about a sun-spot and the corresponding layer rela- 
tively quiet, we shall have large magnetic action and a very small 
resultant electric field over the layers; we must then assume their 
location to be below the layers accessible to spectroscopic investi- 
gations. In this manner we succeed in bringing the electric field 
to a very small fraction of the value given above, but yet not 
small enough to explain the absence of the Stark effect. 

It seems as if the only possible explanation will be the hypothesis 
that the Zeeman effect is due to galvanic currents instead of con- 
vection currents. This hypothesis will be developed in another 


paper. 


XIII. THE MOTION OF ELECTRIC CORPUSCLES IN THE SUPERPOSED 
ELECTRIC AND MAGNETIC FIELDS OVER A SOLAR VORTEX 


There is considerable interest in the study of the motion of 
electric corpuscles over sun-spots because of various applications 
to the theory of prominences and to the theory of the aurora, 
especially the latter, if we suppose with Birkeland and others that 
the aurora is caused by such corpuscles coming from the neighbor- 
hood of sun-spots. 

If we limit our studies to the highest stratum of the sun’s atmos- 
phere, which may be considered as a vacuum, this investigation 
may be made by the methods which I have applied to the theory 
of the aurora. In the following we neglect the action of the 
general magnetic field of the sun. 

Let us first consider the motion along the axis of the whirl. 
This axis being a line of magnetic force, the magnetism has no 
influence on the motion; we have only to take into account the 
electrostatic field. In order to consider the most general case, we 
will take into account the variation of mass with velocity, a 
hypothesis well founded when the velocity is so great that it is 
comparable with the velocity of light. 

We will follow the equations of Schott.? Let e’ be the electric 


*See my memoirs in the Archives des Sciences Physiques et Naturelles, Geneva, 
1907 and 1911-1912. 

2See Electromagnetic Radiation and the Mechanical Reactions Arising from It, 
Cambridge University Press, 1912, Appendix F. 
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charge of the corpuscle in electrostatic units; then the force acting 
on it will, with the same notation as before, be e’E, and we obtain 


8 
)='E= _ Ww ; 


Ve—v bz 


d (= 
dt 


where m is the mass of the corpuscle for zero velocity; c is here the 
velocity of light, and v the velocity of the corpuscle. 
From this equation we find, as Schott has shown, 
om 


SS 8 e’W+C’ ’ 


V e—v* 





C’ being a constant of integration, which can be determined if we 
suppose the velocity given for a given value of z. Let us for 
instance suppose v=o for z=2, and let W, be the value of W for 
Z=2Z. We get 

cm=—e'W.+C’, 


that is, 





= em+e' (W— W.) . 


V —v’? 


Here 





W =ameh[V p2+-2—V p?+2?] 
W.= 2meh[V pi+-z—V p:+z,] 


If we consider the motion over the whirl in general, the mathe- 
matical treatment is especially interesting for the case in which 
the current-lines are circles, as they are when w=o. If we do, 
not then take into account the general magnetic field of the sun 
but consider only the solar vortex, we have a field that is sym- 
metrical around an axis, which affords very important simplifications 
in the mathematical treatment. 

Let us first suppose the velocity of the corpuscle negligible as 
compared with the velocity of light. The corpuscle being in 
motion im vacuo and under the action of a superposed electric 
and magnetic field, we can apply the researches which I have 
published in the Comptes Rendus for March 11, rg1t. 
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The equations of motion are 


dx eT dz dy) 
We Awa a | 


ote eT ., d% dz], , 
map Les As te» 


+¢’E,, 


@z eT, dy 4], . 
moe ha att: 








There the z-axis is the axis of the whirl, H,, H,, and H, are 
the components of the magnetic force parallel to the X-, Y-, and 
Z-axes, respectively, and E,, E,, and E, are the components of the 
electric force. m and e’ are the mass and charge of the corpuscle 
and c=3X10” the velocity of light. ©’ is measured in electro- 
static units. 

Now in our case the components Hr, Hy, and H,, previously 
introduced, have the values 


1 5® 1 6® 
Ae=3y,° smo, Hs=—R 5R? 
and 
x 
H.= 5 Hr, H,=% Hr, 
where 








epv cos 6 
on iP {3 d6dpdg 


is the previously introduced funciion of X and z. 
In the same manner 


x 
E,=,Er,  Ey=ZEr 
and ‘ ‘ 
W W 
Er=— yp» E=—- 3» 


where W is the electric potential, also a function of R and z only. 
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By introducing semipolar co-ordinates, R and ¢, we get, as in 
my paper, 


“mi=C—e'W 
do e’ 

—=wC.——@ 
mR Ut C. 2 


which give further 


@R_1 80 d’z_1 80, 
d@ 28R’ di? 2 8 


(a) +(a) -2 


e’ .\3 
C,-—— © . 
2C Cc € 


Here v is the velocity and C and C’ are two constants of integration. 
If s denotes the arc along the orbit and @ the angle between the 
tangent and the plane through the point of contact and the Z-axis, 
we have 


where 





—— 
- sin 9=R qs 
and as 
Po. 
~ dt 
we get, as in my paper, 
Cane ® 





sin 6= ——__—_——— 
RV 2mC—2me'W 
Hence the orbit must lie within the space defined by 


, 

€ 
C,-- ® 

Cc 





I. 





All 


RV 2mC — 2me’'W 
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The detailed study of these regions is of fundamental importance 
in the discussion of the orbits, as in the mathematical theory of the 
aurora, where corresponding inequalities have been found. 

As has already been said, the equations of motion take another 
form if we suppose the velocity of the corpuscle not to be negligible 
as compared with the velocity of light. 

We have then to use the fundamental equations due to Schott," 
and we get 


“( me “)=5 H.@ q.® +e'E,, 





di\Ve—pdat}) cl ’dt © * dt] 
d(_me dy\ e—,,dx ,, dz), , 
As 4 . € ms dt Hs, \t® Ey, 











d(_m ds\ e[,, dy , dl, , 
HS) Ps Hyg [te es» 


where m is the mass of the corpuscle for zero velocity, and where 
c is the velocity of light. 
From this, as Schott has shown, we get 
mes 


——— =C’—¢2'’W, 
V o—# 





C’ being a constant of integration and W the electric potential. 
If now we introduce instead of the time ¢, a new variable r 


defined by 
dt c 


a V 2—?’? 


the equations of motion can be written 


Gx ef. ds dy), vr 
waa z | B55, Hy 7, tees» 
dy eT ,, dx Ge}. ; 

a | As Ge 4.7, (re &,, 
ely Vg Ole 
ae C | He Hy, | te E; , 








t Loc. cit. 
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where 
/ x / y , 
E,=p Er: Ey =p ER 
and 
_ ‘ sw’ 
Er= = 8R ’ E,= — 8s) 
with 
Pad / 7\2 
Ww’ se ct (C A : ) 
2me'c 


This remarkable transformation is due to Schott. Further, 


find as before 


do e’ 


2 wy pe 
mR, C: Cc 


and, as Schott has shown in his book (published in 1912), 


@R_ i 8P dz_18P 
dr? 28R’ dr 28’ 


dR\? 2\2 
(a) +(@)-?: 
a 
_ eT on o +(e") . 
a mR mc 


drt_y c—y_ = mc 


dc CW’ 


where 





Further, 


which gives 


€ 
ne 
pita ale Ce pet ? 
dt C’—-'W R 
and since 
pa = VCE Wp mic —, 
— C'—e:'W 
we find 
e’ 
dd C C.— -" 





ds V(C—eWy—ma  R_ 


we 
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As R “ie sin 6, we have the conditions 


‘ 
, 


€ 
c <<" 


I= - —,- 
~Y (C—2Wy— mc? R 





=I, 








which define the regions in space to which the orbits are restricted. 

We have here a wide field for research which can be carried 
out by the same methods as those employed for the aurora. 

As I have shown in the Comptes Rendus for March 6, 1911, 
November 25, 1912, and February 1o and 17, 1913, and in Videns- 
kapsselskapets Skrifter, Christiania, 1913, these considerations may 
also be applied to other cosmical problems, e.g., to the theory of 
the solar corona and to cosmogony. 


XIV. CONCLUDING REMARKS. IMPORTANCE OF SIMULTANEOUS 
INVESTIGATIONS OF THE AURORA BOREALIS AND 
SOLAR PHENOMENA 


As we have seen in Section XII, the difficulties in explaining 
the absence of the Stark effect, supposing the Zeeman effect due 
to convection currents, lead to the hypothesis of galvanic currents 
around sun-spots. In the first approximation we can apply the 
theory of currents in plates and try to determine the current lines 
and the lines of magnetic force in space over the plate for com- 
parison with the configuration of Ha flocculi. 

The conclusions based on the hypothesis of convection currents 
have then to be modified. This will be done in another paper, 
which is in preparation. \We shall return to these hypotheses when 
more data regarding sun-spot fields have been collected. 

In this connection a detailed study of the aurora, its situation, 
and its altitude will be of considerable interest, especially if we 
have auroras connected with definite sun-spot phenomena as have 
sometimes been observed, for it is most probable that the electric 
corpuscles that cause the aurora are the same as those thrown 
out from the sun’s atmosphere; and having studied their physical 
properties in the aurora, we can draw corresponding conclusions 
for the phenomena in the sun. This connection between solar 
phenomena and auroras, which has been elaborated by physical 
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methods by Professor Birkeland’ and mathematically by myself,? 
will probably be of importance in future solar work. Furthermore, 
the photographic method of measuring the altitude of the aurora 
by means of two stations connected by telephone which I have 
employed on my expeditions’ in 1910 and 1913 is capable of giving 
very exact information regarding the auroral phenomena, thus 
laying a firm foundation for the conclusions we may draw from them. 
During the coming period of solar activity, special attention will 
be paid to these points of view. 


CHRISTIANIA 
January 1914 


* The Norwegian Aurora Polaris Expedition, 1, sections 1 and 2. 

2 Archives des Sciences Physiques et Naturelles, Geneva, 1907 and 1911-1912. 

3 “Bericht iiber eine Expedition nach Bossekop zwecks photographischer Auf- 
nahmen und Héhenmessungen von Nordlichtern,” Videnskapsselskapets Skrifter. 
I. Math.-Naturw. Klasse, 1911, No. 17, Christiania; and ‘‘On an-Auroral Expedition 
to Bossekop in the Spring of 1913,”’ Astrophysical Journal, 38, 311, 1913. 
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The Making of 
Modern Germany 


By FERDINAND SCHEVILL 


Professor of Modern European History in The Uni- 
versity of Chicago 


Beginning his narration at the time when the close 
of the Thirty Years’ War saw the elder Germany riven 
into small fragments, Professor Schevill shows the 
building process which has resulted in the strong 
central empire of today. The volume is the product 
of ripe scholarship, is historically correct, and being 
written in concise and clear style presents a most 


Early Days in 
Old Oregon 


By KATHARINE B. JUDSON 


These chapters from the story of Old Oregon reveal 
how fascinating the presentation of history can be 
made. In the form of episodes in the history of the 
Northwest arranged in chronological order there is 


| given what is practically a connected series of the 


readily comprehended account of how Germany came | 


to be and what it is today. 
for both class and general reading. Portraits, maps, 
etc. Price $1.25. 


The American City 
By HENRY C. WRIGHT 
Charities, New York 


An examination of the purpose of the city, and the 
reasons which govern its location, how it is gov- 
erned and financed; how it takes care of the prop- 
erty, life, and health of its citizens, and provides for 
the education and instruction of the children. Other 
topics discussed are the problems of transit, housing, 
racial differences, etc., etc. Price 50 cents. 


Sociology 


By JOHN M. GILLETTE 


Professor of Sociology in the University of 
North Dakota 


It is eminently suited | 


Napoleon in His 


main events of Oregon history. The tale of “‘The 
First White Man’s Ship,” for instance, tells about 
the first vessel seen by the Coast Indians and what 
happened. Then follow Captain Cook’s visit, the 
coming of the American traders, the great Fur Com- 
panies, Lewis and Clarke’s journey, etc., etc. Price 
$1.00. 


Own Words 


By H. E. LAW and C. L. RHODES 
Translated from the French of Jules Bertaut 


| This compilation of Napoleon’s views on such sub- 
First Deputy Commissioner, Department of Public | 


| tion. 


jects as politics and administration, his opinion on 
love, marriage, and women, and what he said about 
the things he did and the reasons he gave for doing 
them, is made from his published works and his 
public utterances. In these pages he reveals the 
diverse aspects of his complex mentality, his great- 
ness and his failings; the former being perhaps most 
evident in his genius for co-ordination or accom- 
plishment, and the latter in his lack of moral obliga- 
Not only is the work a valuable addition to 
Napoleonic literature, but it is also of great general 
interest. Net $1.00. 


Our American 


-~Wonderlands 


Sociology attempts to explain that great phenom- | 


enon which we term human society as a collectivity. 
It seeks to make society an object of scientific study. 
This outline or primer of the subject has been pre- 
pared to meet the demands of those who have some 
curiosity about the nature of society and sociology 
but have not a great deal of time in which to satisfy 
it. It is eminently suitable also for study clubs and 
secondary schools as a guide. Price 50 cents. 
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A series of little journeys to the famous wonder 
places of the United States. The Yellowstone, the 


Grand Canyon, Niagara, Lake Tahoe, and many 


| other lesser-known places are described with all that 


enthusiasm for nature’s mighty works for which Mr. 
James isfamous. As a guidebook the work is inval- 
uable, and as a beautiful gift book it will find a ready 
sale. A distinctive feature of the book is the illus- 
trations, which have been chosen with great care. 
Crown 8vo. Net $2.00. 
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viiit+163| 2% the author's own words— 


pages. “This short history of Belgium, of course, is not a ‘war book’ and can- 
tloth, | not be placed among those books classed as war literature. By consult- 
$1.00 | ing this history, however, the reader will be able to 
postage | understand as he never did before why the Belgian 
extra, | nation of today took the stand it has taken in the 
—— great war, and preferred honor in place of dishonor, 

—! and struggle for freedom in place of ease.” 
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**A bird’s-eye view of the history of Japan’’ 


“How often one hears a request for a short history of Japan. Pro- [5 re 
fessor Clement has answered this query in a most admirable manner. | pages. 
This small volume is just the book that a traveller or one desiring the | 12mo, 
facts of Japan’s history in the shortest compass and in readable form cloth; 


will be especially grateful for.” —Dr. John L. Dearing, Yokohama. ——_ 

Every student, every library, every teacher should have a copy at hand = 
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FROM THE PREFACE 
“Tn telling the story of this search for the mode by which the earth came into being, we 
have let the incidents that led the inquiry on from one stage to another fall in with the steps of 
the inquiry itself. It is in keeping with the purposes of this series of booklets that the motives 
which set researches going should have their place with the quests that arose from them 
The final story of the birth of the earth will come only after a time, when the vestiges of creation 
have been more keenly discerned and more faithfully rendered than is possible now.” 
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